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ABSTRACT 



Measurements of the ratio of the neutron-proton 

electric form factors, /G^ , were made from elastic 

n p 

electron-deuteron scattering to a precision of approx- 

2 

imately 1 to 2 % for the range of momentum transfers, q , 

2 -2 

of 0,10 < q < 0,8 F , and for electron scattering angles 

of ^5° to 120°. It was found that within experimental 

errors the slope as determined from the ratio G^ /G^ , 

n p 

agrees with the thermal neutron-electron interaction slope 
when relativistic corrections and proper deuteron wave 
functions are applied to the electron-deuteron scattering 
results . 
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I. INTRODUCTION 



A. EXPERIMENTAL OBJECTIVES 

The charge form factor of the neutron, G„ , as deter- 

n 

mined by electron scattering is found to be essentially 

2 2 

zero in the range of momentum transfers, q , of 0.3 < q 

_2 

<_ 100 P . This fact is in apparent conflict with the 

results obtained by the scattering of thermal neutrons 

from atomic electrons. By methods to be described later, 

measurements of the slope of the neutron charge form factor 

2 -9 -2 

at very low momentum transfers of q = 10 F , show that 



dG T 



n 



d(q 2 ) 



2 

q ^+0 



= +0.0193 ± 0.0004 F‘ 



( 1 - 1 ) 



2 

It is very doubtful that a non-zero slope for Gp at q -*-0 

2 -2 n 

can be reconciled with Gg - 0 at q > 0.3 P 

n 

Drickey and Hand [1] in 1962 measured Gg /Gg for 

2 2 —2 ^ P 

values of q , 0.3 < q < 2.2 P to an accuracy of about 2%. 

Their results yielded G„ - 0, but when reanalyzed by Casper 

£j 

n 

and Gross [2] (1966) using recently developed deuteron 
relativistic corrections and new deuteron wave functions, 

the value of Gg was brought into fair agreement with the 

n 

value of (1-1). Thus the entire problem is further compli- 
cated by the extreme sensitivity of Gg to the different 

n 

deuteron models assumed. Hence, if very accurate measure- 

2 

ments of Gg , for q small, can be made by electron-scattering 
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techniques, a comparison with the value (1-1) could be made, 
which might yield an answer to the "correct" deuteron model 
to be chosen. 

It is the purpose of this experiment to (a) demonstrate 

the basic agreement between the thermal neutron-electron 

slope value of G„ and the electron scattering G^ ; and 

n n 

(b) try to determine which deuteron wave function (model), 

best describes the deuteron. 

B. NUCLEAR STRUCTURE FROM ELECTRON SCATTERING 

Since the early 1950's production of beams of high 

energy electrons by linear accelerators (LINACS) has made 

possible the study of the electromagnetic structure of the 

nucleon. It is expected that the dimensions of any particle 

of mass m is of the order of fi/mc, which for a nucleon is 
_ 1 14 

about 2 x 10 cm. Hence, to investigate the structure in 
this size range it is necessary to use a probe having about 
these wave lengths or shorter. Linacs produce electrons 
with just such wave lengths. The 100 MeV NPS Linac thus 
"sees" only the outer structure or surface of the nucleus. 

Electron scattering has several definite advantages 
over other methods of nuclear structure research. Briefly, 
the important reasons that electron scattering is such a 
useful tool are the following: 

(a) The basic interaction between the electron and 
target nucleus is known; the electron interacts with the 
charge-current density and magnetization density of the 
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nucleus. Since the electromagnetic interaction is rela- 
tively weak, of the order of a = 1/137, the fine structure 
constant, measurements on the target nucleus can be made 
without greatly disturbing its structure. 

This is in contrast to the situation with protons, 
neutrons and alpha particles, since they are strongly inter- 
acting particles, and the scattering mechanism cannot be 
clearly separated from structure effects in the target. 

(b) Photons interact electromagnetically with the 
nucleus, but since a real photon has rest mass zero, it can 
transfer only a definite amount of momentum to the struck 
nucleus. This is in contrast to electron scattering where 
the momentum transferred to the nucleus is limited only by 
the energy lost to the nucleus. 

The probability that an incident electron of total 
energy E is scattered elastically at an angle ©• (i.e., a, 
the cross section for elastic scattering at angle 0-) can 
be calculated precisely on the assumption that the target 
nucleus consists of point nucleons. That means the charge 
and magnetization densities are delta functions. The 
actual cross sections differ from the predicted cross 
section if the charge and magnetic moment are spread out, 
and the difference is usually in the direction of smaller 
values of the cross section since the smeared out structure 
reduces the field against which the particles scatter. The 
deviation becomes greater as the parameter, q, the momentum 
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transferred to the struck nucleon Increases. The (three) 
momentum transfer Is defined by 

q = p f - Pj_ d-2) 

where p.. and p^, are the incident and final values of the 
electron momentum. The effect of structure of the nucleons 
is normally expressed by the use of form factors, G, which 
are believed to be functions of the structure and momentum 
transfer only; 

q2/ 2s_ cross section for scattering from physical particles 
' cross section for scattering from point particles 

( 1 - 3 ) 

Note that similar ideas occur in all scattering situations 
in which the structure of the scattering center has an 
influence, e.g., x-ray scattering form factors reflecting 
the electronic structure of atoms. 

C. ELECTRON SCATTERING FROM DEUTERIUM 

The deuteron, as the only bound two-nucleon state, has 
been the source of much information for nuclear physics, 
due initially to the great interest in the two-nucleon 
problem. Originally the two-body problem was thought to 
be of great importance because of its relative mathematical 
simplicity, and because of the hope that the inter-nucleon 
force was additive, so that knowledge of this system's 
solution would allow extension to more complex nuclei. It 
is now generally assumed that the nuclear force is not 
additive (e.g., the nuclear force, unlike the electromagnetic 



force j "saturates") and hence the two-nucleon interaction is 
not necessarily the best source through which understanding 
of the nature of the nuclear force can be gained. 

Since free neutrons are unstable and decay with a half- 
life of 17 minutes to a proton and electron, deuterium is 
a convenient source for neutrons as targets. For this 
experiment deuterium in the form of solid deuterized poly- 
ethylene was used as a target. Assuming that 

deuteron wave functions are known from theory, neutron 
structure is inferred through a comparison of electron- 
proton and electron-deuteron scattering results. Converse^, 
knowing the nucleon form factors from theory makes compari- 
sons of different deuteron wave functions possible. 

There are several possible outcomes to an electron- 
deuteron scattering experiment: 



(1) 


e + D e ’ 


+ D’ 


(elastic scattering) 


(2) 


e + D e ’ 


+ P + N 


(inelastic scattering, break-up) 


(3) 


Radiative 


scattering 




(4) 


Reactions 


with pion 


production . 


(2) 


occurs only 


above 2.2 


MeV excitation energy and is 



not of interest here; (3) can be calculated; (4) is not a 
problem for incident energies smaller than 150 MeV . The 
theory for process (1) was developed initially by Jankus 
[33 (1956) in the first Born approximation, which should 
be adequate because of the small deuteron charge. The 
deuteron is described by nonrelativistic wave functions 
while the electron is treated relativistically . Jankus 's 
result will be presented in Section II, part B3. 
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Attempts at a relativistic deuteron theory were initi- 
ated by Jones and carried forward by several authors, prin- 
cipally Gourdin, Adler and Gross [4, 5, 6, 7, 8, 9, 10, 11]. 
Results of their attempts will be presented in Section II, 
part B5. 

D. THE NEUTRON-ELECTRON INTERACTION 

The neutron-electron (n-e) interaction is the term used 
to describe part of the electromagnetic interaction between 
the neutron and the electron. This interaction may not only 
be a specific interaction between neutron and electron, but 
between any charged particle and the neutron. Furthermore, 
there exists the possibility of some n-e interaction of a 
nonelectromagnetic (weak) nature. For the present case 
however, results that stem from the interaction of an elec- 
tron with the internal electromagnetic structure of the 
neutron will be termed "the neutron-electron interaction. 

Since the neutron has zero total charge there is no 
Coulomb interaction between a neutron and an electron. 

Both particles have magnetic moments and hence there exists 
a spin- dependent magnetic dipole - dipole interaction 
between them and also a velocity-dependent interaction 
between the magnetic moment of the neutron and the magnetic 
field associated with the convection current of an electron 
in motion. These interactions have been well documented 
and are not of interest here [12] . If there are regions of 
non-zero charge-density in the neutron, (i.e., a charge 
separation exists) then any charged particle probing the 



16 



neutron will experience spin and velocity- independent 
electrostatic forces. These are the forces of interest in 
the n-e interaction. 

At least some charge separation in the neutron is 
expected on the basis that the free neutron possesses an 
anomalous magnetic moment. This contribution to the inter- 
action is called the magnetic or Poldy term [13] • The 
remainder of the interaction (due to charge separation) is 
the result of the fact that a neutron can be part of the 
time dissociated into a negative pion and a proton, 

n <--*■ p + + ir" . 

If the neutron spends 80 % of the time as a neutron and 20 % 
of the time as a proton + negative pion, then an electron 
penetrating the 'neutron’ would see an electric field 
strength equal to that produced by a charge of 0 . 2e , e the 
proton charge. The force resulting from the field would 
have a very short range since the tt ~ would screen the pos- 
itive charge at distances greater than the tt ~ p + separation. 

(The separation can be on the order of the pion's Compton 

-13 

wave length ft/m^C ~ 10 cm. ) The force should be 
attractive and can be represented as an extremely narrow 
potential well. The purpose of all the experiments 
described below was to determine the well depth V , essen- 
tially by measuring the scattering amplitude. 

The first of three techniques to measure the n-e inter- 
action was introduced by Havens, Rabi, and Rainwater [1*J] 
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in 19^7. It consisted of scattering low energy neutrons 
(~0,1 eV) from liquid lead and bismuth, noting that the 
total scattering cross section of an atom includes an 
observable term arising from interference between the 
nucleus and the electrons scattering coherently, while at 
higher neutron energies (~10 eV) this interference term 
is almost absent. Results on liquid bismuth yield 

dGg /d(q 2 ) 2 = 0.0225 ± 0.0070 F 2 . 

n q^->0 

The work of Hughes, Harvey, Goldberg, and Stafne [15] 
depends on the fact that the difference in the refractive 
indices for neutrons at an interface between bismuth and 
liquid oxygen (liquid mirror) comes largely from the neutron- 
electron scattering amplitude. The measurement of the angle 
of total external reflection gives the difference between 
the two scattering amplitudes and hence the n-e interaction. 
Here the results are 

dGg /d(q 2 ) 2 = 0.0200 ± 0.0019 F 2 . 

n q^->0 

The most recent and most accurate measurements have 
been by Krohn and Ringo [16] ( 1 9 6 6 ) utilizing a technique 

of Fermi and Marshall [17] (19^7). It depends on the fact 
that in scattering neutrons, the interference between the 
scattering from the nucleus and from the electrons in an 
atom leads to an asymmetric angular distribution of 
scattered neutrons in the center-of-mass system. This 
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occurs because the atomic form factor depends on the scat- 
tering angle. A summary of the Krohn and Ringo results are 
shown in Table I. 



TABLE I 



Gas 


V (eV) 
o 


dG E /d(q 2 ) F 2 
n 


Statistical 

weight 


Argon 


-3780 ± 250 


0.0196 ± 0.0013 


0.13 


Krypton 


-3800 ± 130 


0.0197 ± 0.0007 


0.30 


Xenon 


-3670 ± 90 


0.0190 ± 0.0005 


0.57 


Weighted 
me an 


-3720 ± 90 


0.0193 ± 0.0004 





It is important to note that all measurements agree that 

p 

dG F /d(q ) is not zero or negative. See Appendix A for 
n dG E 

the details of the calculation of n from V . 
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II. THEORY 



A. THE NONRELATIVISTIC DEUTERON 



1 • The Ground State of the Deuteron 

There are three possible states of the two-nucleon 
system, the di-proton (pp), the di-neutron (nn) , and the 
deuteron (np), of which only the deuteron is known to be 
stable. Experimentally measured static properties are: 

(1) Binding Energy W = 2.226 ± 0.003 MeV 

(2) Magnetic Moment = 0.85735 ± 0.00003 nm 

where lnm = 1 nuclear magneton = eh/2m c and 

P 

nip is the proton mass 

(3) Electric Quadrupole Moment, = 2.82 x 10 -2 ^cm 2 
= 28.2 F 2 , where IF = 10 -13 cm. 

(4) Spin Angular Momentum, S = lfi 

(5) Radius r d = 4.3 F (see page 23 for definition of r^) 

(6) Mass M^ = 2.01410 umu (unified mass units) 

= 9.50 F” 1 



Because the neutron has zero charge, the nuclear force 
binding the deuteron cannot be electrical and neither can 
it be gravitational, for this force is far too weak to 
provide a 2.2 MeV binding energy. Thus a nuclear force 
that is strong and attractive is postulated. By ignoring 
some experimentally determined parameters a simple quantum 
mechanical treatment of the deuteron is possible. 

Assume that the nuclear force is central, i.e. the 
interaction potential between the neutron and proton is 
some function V(r), where r is the separation of the two 
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nucleons. This assumption is in conflict with (3) above 
because a non-zero quadrupole moment implies a nonsymmetric 
charge distribution. Nevertheless, assume the central 
potential, and take the ground state of the deuteron to be 
S^. Hence the Schroedinger equation in the center-of-mass 
system is, 

2 

V 2 ^(r) + V(r) ip(r) = Ei/>(r) (2-1) 



where the reduced mass is 



m = 



mm M 

p n „ M 

m + m 2 

P n 



(2-2) 



and m , m are the masses of the proton and neutron. The 
p n 

average of these masses is M, and r is the relative vector 
coordinate between nucleons. E is the energy of the rela- 
tive motion of the two particles. Since the ground state 
is spherically symmetric (S state), ij;(r) depends only on 
r = | r | , so that 



<K?) 



, , \ u(r) 
'P(r) = — — - 



(2-3) 



The quantity u(r) is called the reduced radial wave function. 
Now using (2-3) and (2-2), (2-1) can be rewritten as 



zhl 

M 



u"(r) + V(r)u(r) = Eu(r) 



(2-4) 



where u"(r) = — and E = -W = -2.226 MeV. 

dr 

square-well potential 



For the 
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V(r) = 



-V 



o 



r < r 

— o 



(2-5) 



r > r 



o 



(2-4) becomes 



Case I r < r 



o 



u"(r) + a u(r) =0 (2-6a) 



where 



a = 



M 



2 ( V E> = 



M 



P ( V -W) 

2 ° ( 2-6b ) 



The general solution of (2-6a) 
is 

u(r) = A sin ar + B cos ar 

(2-7) 

Requiring that i|>(r) be finite 
at r ->■ 0 and i^(r) = 0 at r 00 
forces B = 0 in (2-7), hence, 



u(r) = A sin ar 



(2-8) 



Case II r > r 



o 



u" (r) - y u(r) = 0 (2-9a) 



where 



MW 

Y ~J h2 



(2-9b) 



The general solution of 
(2-9a) is 

u(r) = Ce _Yr + De Yr (2-10) 



the boundary condition at <» 
requires D = 0 so that 



u(r) = Ce Yr (2-11) 



Now at r = r Q both and ip ' (the slope) must be continuous, 
thus (2-8) and (2-11) state 

-yr 



A sin ar = Ce 
o 



o 



(2-12) 



and 



-yr 



Aacos ar Q = -yCe 



(2-13) 



The ratio of (2-13) for (2-12) yields 



a cot ar = y = 
o 



MW 

^2 



(2-14) 
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This is a transcendental equation relating the range of the 
potential r Q to its depth V Q . Typical solutions of (2-1*0 
for r Q vs depth are given in Table II. 

TABLE II 



Range r (F) 


Depth of 
Potential (MeV) 


1.0 


120 


1.5 


59 


2.0 


36 


2.5 


25 


oo 


5.83 



The results of this rough calculation are quite good. 

Independent measurements of r Q yield numbers on the order 

of 1 Fermi. If r = 2.0 F is accepted, the size of the 

o 

deuteron can be specified by 1/y = *K3F, which is about 
twice that of the range r Q of the potential. This explains 
why the deuteron is a loosely bound system (its binding 
energy is 1.113 MeV/nucleon, vs. 8.0 MeV/nucleon in the 
average nucleus). 

There is no experimental evidence for any excited state 
of the deuteron. It can be shown by simple theoretical 
arguments that in order to produce a bound excited state 
the potential well depth must be much greater than in the 
ground state. It must be concluded that no bound state 
exists for H / 0 in the deuteron. 
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2 . Review of Low Energy Neutron-Proton Scattering 



Certain parameters pertaining to the deuteron have 
been established by investigation of the unbound n-p state, 
that is through neutron-proton scattering experiments. It 
is clear that these parameters should be predicted by any 
good theory of the deuteron. With this in mind this brief 
review is given. 

The n-p scattering cross section has been examined 
closely, both at high and low energies. The cross section 
depends strongly on the energy of the incident neutrons. 

At low energies, below 10 MeV, the scattering is essentially 
due to neutrons with angular momentum A = 0 ( S-wave scatter- 
ing) . Thus the angular distribution of the scattered 
neutrons is isotropic in the center-of-mass system. Since 
the proton has molecular binding energy of about 0.1 eV, 
only neutrons with incident energy greater than 1.0 eV are 
considered. Hence the neutrons strike essentially free 
protons . 

The Schroedinger equation for the scattering problem is 

V 2 iJj(r) + ^ [e - V(r )1 ip(r) = 0. 

h 2 L J 

For a central potential and for neutrons incident along the 
z-axis the above becomes 

V 2 ^(r,0-) + — [e - V(r)l i|j(re-) = 0 (2-15) 

h 2 L J 



which has the asymptotic solution 



T|>(r ,£-) 



A 



+ f(e-) 



( 2 - 16 ) 



r» -> oo 



ikz 

e 



ikr 
0 

r 



The f(O-) is called the scattering amplitude representing 

ikz 

the scattered neutrons (waves); Ae represents the incom- 
ing plane wave. It is not difficult to show that 



in - <2-n) 

Using separation of variables, \p(r,&) = R(r) Y(U-), (2-15) 
divides into a radial and angular equation. The angular 
part is 

5TH¥ & ( sln(> H) + K(t+1 > Y - 0 < 2 - 18 ’ 

with solutions Y^(O-) = P^(cos ©•) , the Legendre polynomials 
of order A. Since the potential is central these solutions 
are independent of the potential V(r). The radial equation 
is 



" / \ , ,2 tt/\ A ( A+l ) 

u (r) + k - U(r ) — 

r 



u t (r) = 



(2-19) 



where 

,2 2mE tt/ \ 2mV(r) 

k = — 5 - and U(r) = 5 — L 

fi ft 

and R(r) = u(r)/r. Let the potential V(r) be "turned off." 
Then the solutions of (2-19) are again the spherical Bessel 
functions j^(kr). It can be shown that as r -> 00 

j (kr) - sln(kr- * n/ 2 ) „ y off „ 

36 Kr 
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and hence 



u A (r) — v oo - sin(kr - Air/2). (2-20) 

When V is "turned on" but r is still greater than r , the 
solution of ( 2 - 19 ) must still correspond to that of a free 
particle, and hence be the same as (2-20) except for 
possibly a phase factor. Thus 

U l/ r ^~ r sln ( kr ~ ^jr + ^ " v on "- (2-21) 



The general solution to (2-19) for "V on" at r > r Q 

oo 

♦ (r.eo = £ Cj Pj[ ( cos& ) 

t = 0 



is 



( 2 - 22 ) 



■ I 

£=0 



cos©-) . 



But this must be the same as (2-l6). Now using a Rayleigh 
i k z 

expansion on e in (2-16), and writing all quantities as 
complex exponentials, equating (2-16) and (2-22) yields 

„ _ i (2t+l) e 

C £ " k 



and 



f ( ©•) = (2£+l)e ^ sinfi^P^ ( cosO-) . (2-23) 

l 

The square of the scattering amplitude is the differential 
scattering cross section (2-17), hence 



O r 



■/ 



do(e-) 

dq 



sinG- dO- 



jhr V 

k 2 ^ 



( 2£+l) sin^6 



(2-24) 
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Note that S-wave scattering (£ = 0) is predominant for 
energies less than 10 MeV. By taking the range of nucleon 
force to be 2F and assuming a 0.5 MeV incident neutron 
(in center-of-mass frame or 1 MeV in the lab) it can be 
seen by forming the ratio of the coefficients in (2-22) 
that 



B 

B 




(kr Q ) 



2 



This states that only about 9% of the scattering is due to 
neutrons with l = 1. 

In (2-24) let l = 0, then 

h_ 0 /sin 6 \2 

°T = jL. sin 2 « o = ■ <2-25) 

As k -> 0, sin 6 q approaches S q . In order that a,p remain 
finite <5 0 /k must approach some constant, called -a by 
convention. Thus 

a T = 4ua 2 (2-2 6) 



where a is called the scattering length. Consideration of 
the solution 

6 

U Q (r) - sin(kr + 6 q ) = sink(r+-^-) k _ vQ '“k(r-a) 

demonstrates the fact that the asymptotic solution becomes 
a straight line, and that the scattering length, a, is the 
intercept on the r-axis and is obtained by extrapolating 
the radial wave function inside the potential, well beyond 
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the range of force r Q . Figure 1 illustrates the signifi- 
cance of the sign of the scattering length. The scattering 
length a is positive if the state is a bound state. 





(a) Figure 1 (b) 

Positive scattering length Negative scattering length 
indicates a bound state indicates an unbound state 

The energy dependence of the scattering cross section 
can be expressed in terms of the scattering length a, and 
another parameter which has the dimensions of length and 
is known as the effective range. Denote this parameter 
by p, where p is of the same magnitude as the range of 
nuclear force r . Effective range theory predicts the 
phase shifts 6^ as a function of energy. Consider S-wave 
scattering only. The wave equations for two energies 
E = h^k^/2m and E = 0 are 
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Case I V / 0 

u" + (k 2 - ^~)u = 0 (2-27a) 

fi 

u" - = 0 (2-27b) 

° ft 2 

multiply (27a) by u Q , (27b) by 
u and subtract, getting 



^ (uu i - u o u ’) * k2uu o (2 ‘ 28) 



Subtract and integrate the above 



Case II V = 0 
v" + k 2 v = 0 (2-29a) 

vj = 0 (2-29b) 

multiply (29a) by v 0 , (29b) 
by v and subtract, getting 

^( V v;-VoV')=k 2 v o v (2-30) 
equations. Then 



(uu 1 - u u' - vv' + v v ' ) 
o o o o 



p=oo 



r=0 



= k 



/ (u %- 



vv o )dr. (2-31) 



Because u = v, u = v for r > r , the expression on the 
J o o o’ K 

left hand side contributes nothing at the upper limit. 
Furthermore u(0) = u Q (0) = 0, hence (2-31) is 



<VV o - 



v c v') 



= k‘ 



r=0 



oc 

1 



(uu -vv )dr. 
o o 



( 2 - 32 ) 



In analogy with (2-20) the solutions u and v can be written 

sin(kr+6 ) 

“Cr)— v( r) = C sin(kr+« 0 ) = sln <i o 

where C is a normalization constant. Now (2-32) becomes 

00 

- ^ - k cot 6 q = k 2 / (uu o -vv Q )dr 

and in the limit of zero energy neutrons u - u q , v * v Q so 
that 

-1 k 2 

k cot 6 = "a + P 2 
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where 



P 5 2 



oo 

/ 

o 



( v 2 - u 2 )dr 
o o 



is the effective range. 

Hence the cross section is given by 

4tt 



°T 



k 2 + ( i_ PlL 2 
k C a 2 



(2-33) 



-) 



2 2 

Besides being a function of k = 2mE/ft , the cross section 
is expressed in terms of the two parameters, effective range 
p and scattering length a. The cross section is independent 
of the . form and shape of the potential. Both p and a can 
be determined by measuring the cross section at different 
low energies of the incident neutron. 

The theoretical estimate of the value of the scattering 
cross section for* i = 0 is far below the experimentally 
observed value. In fact for a square well of reasonable 
width and depth 

- 4.5 x 10 -2 ^ cm 2 

while 

0 . , , = 20.4 x 10" 2iJ cm 2 , 

experimental 

This conflict was removed when Wigner realized that the 
theory for is based on a potential whose "strength" is 
calculated from the deuteron, where the spin is 1, i.e. 
neutron and proton are parallel, while scattering experi- 
ments are done using randomly oriented nuclei. The 
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scattering experiments then represent an average over 
parallel (triplet) and anti-parallel (singlet) states. The 
weight factor for the singlet and triplet cases are 1/^1 and 
3/^ respectively, hence the total cross section for scatter- 
ing is 

^ = i a t + } a s 

where a ^ and c g are the triplet and singlet cross sections. 

Using (2-33) and working backwards assuming = 20. *1 barns, 

1/Y{- = ^.3 P, and r^. - r g = 2 P it can be determined that 

W = 0.066 MeV with 1/y = 25 F, where the scattering length 

s s 

a is negative. Hence the singlet state is not a bound 
deuteron state. 

3 . Pheonomenological Two-Nucleon Potentials 

On the evidence presented so far it can be concluded 
that: (1) The deuteron is predominantly in a state 

since the solutions to the Schroedinger equation for a 
central potential are quite good; (2) since ^ + y^, 

and ^ 0 , there is reason to guess that noncentral forces 
are acting (the inequality of magnetic moments in a complex 
problem and could well have other interpretations); (3) the 
forces between neutron and proton are different in the 
singlet and triplet states, implying that nuclear forces 
are spin-dependent. There is a point of confusion on (3). 
Spin-dependent does not imply that this contribution is not 
central in nature. It means that the usual central scalar 
potential is different according to the relative orientation 
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of spins, a situation that has no direct analogy either in 
electromagnetism or atomic physics. 

A realistic deuteron potential must incorporate these 
ideas and more. The potential must be symmetric in the 
coordinates of the two nucleons and invariant under rota- 
tions and reflection of their space coordinates. This 
means the potential is a scalar. 

Some of the more detailed features of the nuclear 
potential can be directly related to the spin and isospin 
dependence of the nucleon interaction. An important 
contribution for the nuclear symmetry potential comes from 
the exchange character of the nuclear forces. Thus a 
Serber exchange potential, P r , for the central interaction 
would be 

V c = | (1 + P r ) V(r) or 

V Q - | [1 - | (1 + t 1 -t 2 )(1 +a 1 *a 2 )]V(r) (2-3*0 

-)■ -)■ — V 

where t^, and a^, o^ are the isospin and Pauli spin 
matrices respectively. 

The occurrence of a rather strong spin-orbit force in 
the nuclear interaction gives rise to a spin-orbit coupling 
in the two body force; 

V LS = h V LS^ ^ r l " r 2' ) X ^ P 1 “ p 2^ ’ ^ S 1 + s 2^ ' (2-35) 

This V^g contribution is based on Wigner's idea that 
velocity-dependent forces are acceptable when they depend 
only on the relative momentum p of the two nucleons. 



32 



All these facts taken together yield a general potential 
V = V c (r) + V a (r) o 1 'o 2 + V T ( r ) S 12 + V LS ( r ) L • S (2-36) 

where a^*a 2 re P resen ts the spin-spin Interaction and S^ 2 is 
the tensor operator 

S 12 = 3(a 1 -r ) (a 2 *r) - (a^Cg). (2-37) 

The form of the tensor operator can be inferred by analogy 
with the classical problem of the potential energy of an 
electric dipole d^ in the field of another dipole d 2 at a 
distance r.* Here 

V dd ~ r ~ 3 [3(d 1 *?)(d 2 -?) - d 1 *d 2 ]. (2-38) 

The r comes from the fact that the electric charges 
interact with a potential proportional to r - '*". But the 
r-dependence of "mesic charges" which are thought to be the 

_-5 

source of nuclear fields is unknown. Hence the r J depen- 
dence must be replaced by an arbitrary function of r, V,p(r). 

The structure of the nuclear force is closely related 
to the properties and interactions of the whole family of 
strongly interacting particles. It has not been possible 
so far to derive the forces between these particles, or 
their masses, from simple assumptions regarding the basic 
structure of the strong interaction. Nevertheless, certain 
relations between the interactions and masses of the 



*See Appendix B for properties of the tensor operator 

S 12 * 
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strongly Interacting particles can be established. Of 
particular significance is the relationship between the 
nuclear force at large distances and the pion-nucleon 
interaction . 

The force at large distances can be described in terms 
of an exchange of mesons, in a similar manner as the elec- 
tromagnetic interaction can be analyzed in terms of photon 
exchange. The special role of pion exchange is due to the 
smallness of the pion mass as compared with that of other 
strongly interacting particles. The interaction associated 
with the exchange of a particle of mass is limited to 
a range of the order of the compton wave length A c =fi/m^c. 

In fact, the intermediate states involved in such an 

2 

exchange have an energy of at least m^c and so are limited 

2 

to a duration of order tfi/m c . During this time, the 
emitted particle cannot travel farther than h/m^c and thus, 
the interaction is expected to decrease strongly for 
distances greater than A . If n particles are exchanged 
simultaneously, the corresponding intermediate energies 
are n x m^c and the range is A /n. The interaction at the 
largest distances is thus determined by the exchange of 
single mesons. 

The asymptotic form of the one-pion exchange potential 
(OPEP) is uniquely given by the mass and symmetry properties 
of the pion (J 71 = 0~, T = 1) 



3 ^ 



/ 






V, 



OPEP 




+ ~2‘)] i T" (2-39) 



where x = m^ c/h = 0.70 F - 1 , f 2 /hc = 0.081 ± 0.002. 



The strength of the potential is determined by the coupling 
constant f for the process N -*• N + it. The meson field is 



Equation (2-39) reduces to (2-36) for the deuteron. 

Early definitive work on the phenomenological nucleon- 

nucleon potential in the OPEP approximation was carried 

out by Glendenning and Kramer [18] . They attempted to 

construct triplet (spin) - even (parity) potentials that 

are asymptotic to the OPEP and are modified in the inner 

regions so that the deuteron properties are obtained. The 

results of an 8-parameter fit to an equation similar to 

(2-39) are found in Table III, entries 1-9. Coefficients 

in their potential are varied such that the empirical values 

of = 28.2 F 2 and a^ = 5-^ ± 0.02 F are reproduced. A 

hard core radius r was chosen so that the potential has a 

c 

bound state at W = 2.226 MeV. (The hard core is discussed on 
page 3 8). 

Many other authors have contributed to the data [19, 20 , 
21] . One of the most useful has been the Hamada-Johnston (HJ) 
[22] potential given by 



7T — 

pseudoscalar since J = 0 , and isovector since T = 1. 




(2-110) 
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TABLE III 



Potential 

Identification 


r (f) 
c 


a^F- 1 ) 


P(F) 


Q d (F 2 ) 




X 


GK1 


0 .4815 


5.376 


1.717 


28.79 


7.415 


0.02714 


2 


0.2466 


5.368 


1.709 


28.21 


6.281 


0.02693 


3 


0.4747 


5 . 456 


1.812 


28.17 


5.970 


0.02595 


4 


0 .4425 


5.364 


1.703 


28.22 


7.103 


0.02677 


5 


0.3632 


5.366 


1.705 


28.02 


6.710 


0.02663 


6 


0 . 3924 


5.384 


1.726 


28.12 


5.957 


0.02676 


7 


0.5369 


5.477 


1.836 


28.14 


6.029 


0.02564 


8* 


0.5007 


5.413 


1.760 


28.18 


5.622 


0.02654 


9* 


0.4329 


5.373 


1.715 


28.29 


7.425 


0.02670 


HJ 


0 .485 


- 


1.77 


28.50 


6.97 


0.02656 


FL 


0.735 


- 


- 


- 


4.60 


- 


Experiment 


- 


5.40 + 


- 


28.20 ++ 


- 


- 



X is the asymptotic D state/S state ratio 

* GK8 and 9 include the spin-orbit term in (2-3*0 while GK1-7 
are without it . 

t N. K. Glendenning and G. Kramer, Phys. Rev. 126 , 2159, 1962. 
++ G. P. Auf f ray , Phys. Rev. Letters 6 _, 120, 1961. 
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where 



L 12 = L 2 -|-[(a 1 -L) (a 2 *L) + (a 2 *L) (a^L)] 



= (6 lj + a 1 -a 2 ) L 2 - (L-S) 2 



and contains central V c , tensor V T , spin-orbit V LS , and 
second-order spin-orbit V LL components. It is from this 
potential that the later used Partovi wave functions have 
been derived. The radial functions are restricted by the 
condition that at large distance the central and tensor 
potentials should be described by the OPEP. The radial 
forms are : 



V c = v q (t 1 *t 2 ) (a 1 -a 2 ) Y(x) [1 + a c Y(x) + b c Y (x)] 
V T = v o < ^l’^2^l’^2' ) Z(x) [1 + a T Y(x) +b,pY 2 (x)] 



V LS S LS v o Y ^ x ^ 1+b LS 



V LL = «LL v o ^ [l+a LL Y(x) + b LL Y 2 (x)l 



X 



1 f 2 ? 

v = i me =3.65 MeV 
o 3 he tt 

me 

it r v/ N e 

x = — — r = - — 1 , Y(x) = - 



-x 



h 



1. 43F 



X 



Z(x) = 1 + ^ . 

x 



In addition the potential has been assumed to have a com- 
ponent giving rise to an infinite repulsion at radius 

x c = 0.343 (r = 0.485 P) . 
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The values of the quantities V , V,p , etc. have been 
calculated for four different types of potentials. The 
results are shown in Table IV. 



TABLE IV 







V c (MeV) 


V T (MeV) 


V Lg (MeV) 


V LL (MeV) 


Singlet 


even 


-1460 


- 


- 


-42 


Triplet 


even 


- 207 


-642 


34 


668 


Singlet 


odd 


2371 


- 


- 


-6688 


Triplet 


odd 


- 23 


173 


-1570 


-1087 



Values of HJ potential at r = 0.485 F. 



The "hard core" mentioned above is a phenomena that 
cannot be investigated through low energy nucleon-nucleon 
scattering. In the neighborhood near and above 200 MeV 
repulsive components in the nucleon-nucleon interaction are 
seen. One way of introducing a repulsive interaction (which 
might account for nucleon saturation) is in terms of a very 
strong repulsive potential of very short range - "the hard 
core." The infinite short ranged repulsion is only one, 
rather extreme, way of accounting for the observed effects. 
Finite repulsive potentials, "soft cores," as well as inter- 
actions depending explicitly on relative velocity can account 
for some experimental results. 
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A fairly recent idea by Feshbach and Lomon (FL) [23] in 
1967 > is to replace the interior region by appropriate 
boundary conditions at some finite radius. As long as the 
interactions responsible for these boundary conditions 
involve energies much greater than the bombarding energy, 
it may be expected that the boundary conditions will be 
approximately energy independent. The PL model also employs 
the (TPEP) two pion exchange potential. 

4 . Solutions of the Non-Central Force Deuteron Problem 
Solutions to the non-central force deuteron problem 
are obtained by substituting the potential Vq PE p into the 
Schroedinger equation and solving for the wave functions. 
Angular momentum and parity arguments can determine the 
deuteron ground state angular wave function. 

The total angular momentum of the deuteron is J = 1 
and it is in a definite parity state. The value J = 1 can 
be obtained from different combinations of the orbital 
angular momentum l, and the spin s such as: 



H = 0 
1=1 
1=2 



s = 
s = 
s = 
s = 



1 

0 

1 

1 



3 

1 

3 

3 



S 

P 

P 

D 



1 

1 

1 

1 



and no others. If the deuteron is a mixture of these it 
3 3 13 

can only be S^ + or P^ + P^ since the combination 

must have a definite parity. Since an almost spherically 
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symmetric ground state exists, the obvious choice is the 
S state combination. It is not difficult to show that the 



1 ? 

odd parity states ( yield the wrong values for 



Q d and y d< 



a , M 

The angular momentum wave function 1A for a two 

particle system in the orbital angular momentum state l and 
spin state S, forming the total angular momentum J, with 
Z-component M is given by 

- 2, <1Sm,,M s |jM> x s (2-J41) 



H 



M 

&SJ 



M 0 M 

io S 

where the are the spherical harmonics, x s are the spin 

wave functions and <& Srru M |JM> are the Clebsch-Gordon 

X/ s 

coefficients. Two spin ^ particles have the following spin 
states : 



Triplet (s = 1) 

Xi = K hhhh\ ll>a(l)a(2) = a(l)a(2) 

X° = <l &rhh\ 10>a(l)3(2) + <hyk~h\ 10>a(2)B(l) 




( 2-^2a) 



(2-*J2b) 



X 




Singlet (S = 



X 



0 

0 



X 



0 

0 



<hh~h~h\ 1-1>B(1)3(2) = B ( 1 ) 3 ( 2 ) (2-^2c) 

0) 

00>a(l)B(2) + <¥s“¥s| 00>B(l)a(2) 

— fa(l)B(2) - $(1 )a (2)1 (2-^3) 

/2 I J 
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where a = (q) and 6 = (^) are "spin up" (proton) and "spin 
down" (neutron) states of the nucleon. The arguments of 
the a's and 3’s stand for particle number one or two. The 
angular momentum wave functions "Ij. for the two states 

under consideration are 

J = 1, M = 1 j & = 0 j s = 1 

= <010l|ll> Y°a(l)a(2) = Y°a(l)a(2) (2-44) 




1 

Oil 



J = 1 j M = 1, 1=2, s = 1 



1 

211 



^ 1-M M 

= 2 <21(l-M s )M s |ll> Y 2 s Xl s 



M 



= yf [ Y 2 3(1)3(2) Y 2 -^[a(l)B(2) + 3(l)a(2)] 



+ Y° a(l)a( 2) . 



( 2 - 45 ) 



is 



The Schroedinger equation in the center-of-mass system 



where 



v 2 + v c + Vl2 



4> = EcJ> 



(2-46) 



4 > = d> + = -“ -( - g l 

T Y s t D r 



y l , w(r) 1 / 

oil + -7-f 



1 

211 



(2-47) 



is the ground state deuteron wave function.' For ease the 
linear spin-orbit and quadratic spin-orbit parts of the 
potential (2-40) have been ignored. Rearranging (2-46) 
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tf v2 + v o - E 



on + 



oil 



W 



+ V T (? «%d * r 



;[^oil 



211 



= 0. (2-48) 



In writing (2-48) the following properties of S.^ have been 



used : 






1 
211 



(2-49a) 



and 
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U 1 . 



^ Toil 2 / 211 ' 



(2-49b) 



Note that 



„2 _ 1 9 , 2 9 s L 

r 2 9r (r " ^2 



( 2 - 50 ) 



where 



2 _ -h £ 



L = 



sin 2 ©- 



• n 9 / n 9 \ , 9 

sin 9- (sin 0- g^) + — 2 



(2-51) 



and 



? m m 

ITY £ (e-,(j)) = l(SL+l)fi Y £ (e-,<|>). 



(2-52) 



Hence 



V 2 ( — 'll ■*■) 

^ r /Oil ' 



1 d 2 u 

r dr 2 /( Oil 



(2-53) 



and 



’ 2 <?y 2 n> 



/ 1 d w 1 2(2+1) 



r , 2 r 
, dr 



w 



' 211 ' 



(2-54) 



Multiply (2-48) by oil* inte S rate over the solid angle, 

y M 
£SJ 

to get 



u"(r) + (V c - E)u(r) = -/8 V T w(r). (2-55) 

Similarly, multiply (2-48) by 2 / and integrate to get 

"fnr [ w " (r> - 6 + [ v c - 2V T - E ] w(r> 

= - /8 V T u(r). (2-56) 



These coupled differential equations were first obtained by 
Rarita and Schwinger [24] . 

A method for solving these equations has been described 
by Hamada [21] , who states that (2-55) and (2-56) have two 
sets of linearly independent solutions. Thus choose (u^,w^) 
and (u^jW^) with the asymptotic requirements that the first 
solution behave as an S wave function and the other as a D 
wave function: 



and 




-ar 



e -ar [l + 3/ar + 3/(ar)^j 



(2-57a) 



(2-57b) 



Starting with the behavior at large r, integrate (2-48) 
inward to get the two solutions for arbitrary r. Then u(r) 
and w(r) are obtained as a linear combination of these two 



solutions : 



( 2 - 58 ) 



( 




1 

(l+\ 2 ) h 



/ U 1 ( r) 
■\ W l( r ) 



+ A 



u 2 (r)\ 
w 2 (r)/_* 



The requirement that u and w in (2-58) vanish at the core 



radius r gives an equation for r which is solved by an 
c c 

interpolation procedure. The asymptotic D state/S state 
ratio X is then determined from 



X = - 



u l< r c> 

u 2 (r c ) 



Vrp 

W ' 



(2-59) 



The solutions (2-58) are found in practice by numerical 
methods. Now having the wave functions in hand the magnetic 
moment and quadrupole moment can be predicted by taking the 
expectation values 









<t> dr, 



u = — i + ks +es 
z 2 z to p pz to n nz 



and 



< Q d > 




Q d 4> dr. 





2 



r 



(2-60) 

( 2 - 61 ) 

( 2 - 62 ) 

(2-63) 



Typical results are listed in Table III for various poten- 
tials. (See Appendix C for example calculations of < V*„ > 

CJ 

and <Q d > • ) 

B. THE SCATTERING OF ELECTRONS FROM PROTONS AND DEUTERONS 
Traditionally electron scattering cross sections have 
been handled through the techniques of first order 






time-dependent perturbation theory, the main result of which 
is the Fermi Golden Rule No. 2. This rule relates the 
transition probability per unit time, r, to the square of 
the matrix element for the transition. Specifically, 

r - if l<f|Hj|i>| 2 If (2-6H) 

where dN/dE is the density of states like | f > around the 
energy E, and is the interaction Hamiltonian. The cross 
section is then defined as 

p 

0 incident particle flux' (2-65) 

Quantum electrodynamics (QED), which is the interaction of 
the Dirac spinor field with the electromagnetic field, pro- 
vides the mechanism for computing the cross section via 
S-matrix theory. In general terms the S-matrix is a unitary 
operator that transforms the initial state |i> into the 
final state | f > , i.e., S|i> = |f>, where 

oo 

S fi = I H7 (-D n AV--pV[H I (x 1 )...Hi(x n )] 
n "° ( 2 - 66 ) 
and t is the chronological operator that orders the inter- 
action with respect to causality. The expression (2-66) is 
integrated over all four space and hence is not physically 
measurable. If is normalized to unit space-time volume 



then it can be written 



ij 

where 6 is a four dimensional delta function expressing the 
conservation of four momentum in the transition from | i> to 
| f > . The cross section for a transition from a two-particle 
initial state |i> into some group of final states is in the 
center-of-mass 




1 




E(p 1 ) E(p 2 ) 




final final 

spins momenta (2-68) 



( Integration over suitable variables yields the desired dif- 
ferential cross section do/dft. 

At the heart of (2-68) is the matrix element |M^|. It 
is here the character of the interaction is found while the 
other factors represent kinematical quantities. The form 
of can be surmised by examination of the Feynman diagrams 

for the process under consideration. 

Much of QED is conveniently described by Feynman diagrams. 
Here the diagrams will be mainly used to supplement the 
understanding of the electromagnetic interaction, but their 
use is certainly not limited to this role. If the diagram 
can be drawn, then there exist specific rules for construct- 
ing the S-matrix therefrom. 

Electron scattering is graphically represented by 
Figure 2. Both graphs are equivalent, but (a) is most often 
seen since it is just a little more suggestive of the 
scattering process in the center-of-mass frame. 



(a) 



(b) 



Figure 2 

The time axis is directed upwards unless otherwise noted. 

The horizontal axis can represent a space or momentum inter- 
val, the distinction between the two being of no importance 
here. The arrows on the lines point in the direction of 
increasing time for a particle. If the arrows were reversed 
the particle then represented would be the antiparticle of 
the original one. Figure 3 sometimes is said to represent 
particles in the ’’scattering channel" while the "crossed or 
annihilation channel" is represented by the same diagram 
rotated by tt/2 . 

rotated ccw 
by 

TT 
Z 



scattering 

channel 

Figure 3 





annihilation 

channel 



^7 



Lines beginning and ending within a graph represent "virtual 



particles." Table V below lists some of the more important 
graph elements seen. 

TABLE V 



Graph Factor in the 

Element S-Matrix 



Interaction 





annihilation (creation) 
of an electron 




annihilation (creation) 
of an antielectron 




e u (g) 

v^r 



annihilation (creation) 
of a photon 



iYp - m 

2 2 
p + m 



. 2 , 2 

k + m 



electron progator 



photon progator 



meson progator 




Yj/Cp-P'-k) 



Y P O y <$(p-p’-k) 



photon-electron vertex 



photon-nucleon vertex. 
Bubble or cross hatched 
circle represents ignor- 
ance of the true inter- 
action; mathematically 
this ignorance is repre- 
sented by form factors 



1 . Elastic Electron-Proton Scattering from Point Protons 



To calculate the elastic electron-proton (e-p) scat- 
tering cross section from point (Dirac) protons in the first 
Born approximation it is assumed that (1) the electron is a 
physical point, and (2) the one photon exchange process des- 
cribes elastic scattering. The graph of Figure 2 is then valid. 




From the figure define the kinematical quantities p^jP^, 
Pi>Pf as the four-momentum of the incident and final 
electron and proton respectively ; q is then the four- 
momentum transfer 

q 2 = ( P ± - P f ) 2 = (Pj_ - P f ) 2 . (2-69) 

Choose the convention that q <0 for electron scattering. 

By looking at Figure 4 it can be seen that 



*The metric used is q = q - q . In the lab frame the 
recoil energy given the struck u particle is small^with respect 
to the total incident electron energy thus ^ -q2 < o. The 
four vector q2 is then said to be spacelike. In the annih- 
ilation channel q2 >_ i|M2_^> 0 and is timelike. Another 
metric commonly used is q2=q2_q2 The different metrics 
yield cross sections and form factors with some sign dif- 
ferences. This possibly ambiguity is being remedied by use 
of the variable t for ±q2 in many newer papers. 



( 2 - 70 ) 



s n - M n = j„ 4 jU 

q 

which represents the interaction of the electron four- 
current j with the four-current of the proton J y via 
virtual photon progator ^-2. 

Equation (2-70) is the equivalent of taking the inter- 
action Hamiltonian to be 

H I (x) = j y (x) A y (x) (2-71) 

where A (x) is the electromagnetic field of the proton, 
which in turn depends on the proton current J y (x). Regard- 
less of the point of view the cross section (2-68) becomes 



0 - 6 (Pf-Pi + P f- P i>IJ u ^2 



rk I 2 



( 2 - 72 ) 



q 



QED predicts the form of the current for the electron, 

j y u(p f )Y y u(p j _) (2-73) 

and by analogy for the proton 



J y -U(P f )y y U(P i ) 



(2-7*0 



where u,U and u,U are the incident and final electron and 
proton spinors, respectively; and Y is a *) x ^ matrix. Sub- 
stitution of the currents into (2-72) and carrying out the 
trace calculation of |M fi | followed by a sum over final 
spins and an integration over final momenta will finally 
yield the result 



50 




Dirac 

Proton 



where 




2 — 
cos 2 

. 4 0 - 
sm — 




2 E 
MC 2 



sin 



2 e- 



i-l/ 2 2 „ 

h q ^ 2 0- 

!^ tan 



1 - 



(2-75) 



/ei \ 2 cos 2 I 
\ 2E ) sin 11 | 



Mott 



and 





cross section 



recoil factor 






s 



( 2 - 76 ) 



(2-77) 



are the results of the kinematical quantities in (2-68). 
The square of the invariant matrix element l^fjj yields for 
a Dirac proton 




1 - 



-•2 2 

iLR 



2m 



. 2 0- 
2 tan j 



( 2 - 78 ) 



In obtaining (2-75) the approximation m g -*■ 0 since m g « E 
has been used. Equation (2-75) is correct for scattering 
from a point proton with no anomalous magnetic moment. 

But the proton is known to have an anomalous moment so 
the structure of M^ must be changed to reflect this fact . 

In addition, the possibility of nucleon structure should 
be taken into account in (the electron remains a point). 

Figure 4 is redrawn to schematically illustrate the change 
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Figure 5 



The bubble at the photon-nucleon vertex represents the non- 
locality of the photon-nucleon interaction and the anomalous 
moment of the proton. The matrix element is now 




y/here 



<J y ~U(P f )O y U(P ± ) 



(2-79) 



and 



0 P = F n (q") Y M + F 0 (q")c^q. 



.2w,P 






( 2 - 80 ) 



The quantities F-^ and F^ are undetermined real functions of 
q^, while a y = — {y y ,Y V )- Foldy [25] has shown that while 
(2-80) is not unique it has the proper general form. 

With the above changes (2-75) becomes 



do 

dp 






X £ X 



Mott 



1 2(F 1 +k p P 2 )2 ^ f +K P P I1 

( 2 - 81 ) 



^m c‘ 
P 



where the final bracket is the invariant matrix element 



|M fl | , and where = 1.79 is the protons anomalous (Pauli) 
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moment. Equation (2-8l) is known as the Rosenbluth cross 
section. The factor F^(q ) is found to represent the 
proton's charge distribution and its normal magnetic moment, 
while P 2 (q ) represents the anomalous part of the magnetic 
moment. F^ and F 2 are called the Dirac and Pauli form 
factors respectively. Yennie [26], and later Hand, Miller 
and Wilson [271 pointed out that F^ and F^ should not be 
considered to be electric or magnetic form factors. More 
properly certain linear combinations of F^ and F^ defined 
to be 



°E = P 1 + 



r-2 2 

■ft q „ 

2.2 K p F 2‘ 



4m 



(2-82) 



and 



°M ■ P 1 + VV 



p 

q -»■ 0 



1 + k = V 
P P 



(2-83) 



are termed electric and magnetic form factors. The major 

advantage of (2-82) and (2-83) is that in the Rosenbluth 

2 

cross section (2-8l) the cross term (F^ + K^Fg) is removed. 



da = / da\ 



dQ 



dft J 



Mott 



l 



7 



g 2 _ ft q ' 

E ^ 2 c 2 _ 2 * q 2 

2 2 „ 2 2 2 ta 2 m 



1 



_ £ q 

4m 2 c 2 



4m c‘ 



(2-84) 



In both the cross section formulas the form of the 

2 2 2 0 * 

quantity in braces is A(q ) + B(q ) tan It can be 

shown that this is the most general form for the one photon 
exchange process, regardless of the structure of the 
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electron and the proton. The Interpretation of A and B In 
terms of G £ and G m or F^ and Fg depends of course on the 
structure assumed. 

2 . Nucleon Form Factors 

In the text above no real attempt has been made to 
coherently explain the form factor concept. So far they 
have been characterized mainly as correction factors that 
when multiplied by the scattering cross section for point 
(structureless) particles yield the cross section for 
scattering from physical particles. Thus, as introduced, 
these dimensionless quantities are meaningless unless some 
predictions can be made from a basic theory of matter about 
the structure and behavior of the form factors. 

While the simple idea of the form factor representing 
some physical extension in space can be illustrated by non- 
relativistic examples, a more comprehensive discussion of 
their meaning needs the help of the so-called dispersion 
relations. Both viewpoints are investigated below. 

In Section II, part A2, the scattering amplitude f(0-) 
was found in a partial wave solution of the Schroedinger 
equation for the scattering of neutrons from protons. Of 
Interest here is the solution of the Schroedinger equation 
for the scattering of electrons from two potentials V(r), 
one with a point charge source, the other with a smeared 
out charge distribution. 

With all assumptions as in Section II part A2 , the 
Schroedinger equation for scattering is again 
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(V 2 + k 2 ) iJj ( r ) = 



(2-85) 



2m 

* 2 



V(r)ip(r) 



This differential Schroedinger equation is solved formally 

for ^ making (2-85) into an integral equation. The inverse 

2 2 

of the linear differential operator (V + k ) is an integral, 



hence 






(r) = / G(r ,r ' ) ~ V(r' )^(r' )dr' 

J ft 2 



( 2 - 86 ) 



where r' is a dummy variable of integration and G(r,r') is 

a function characteristic of the form of the differential 

operator. G(r,r') is known as a Green Function. The 

general solution to (2-85) is the sum of the solution to 

the homogenous problem (in this case it would be the 

unscattered plane wave e^ r since ^ V(r) = 0) and some 

h 2 

particular solution to the inhomogenous problem which is 
given by (2-86). Hence 



,i, / s ik*r . 
Y(r) = e + 



Ja(r,r' 



2m 

b 2 



V(r' )^(r ' )dr ' . 



(2-87) 



This is called the integral equation for scattering, and 
the solution is found by iteration. The first iterated form 
is called the Born or plane-wave approximation 



* = e lk ' r + 



2m f ' , 

?r (r - 



, \ ik • r , . 
r ' ) e dr' 



( 2 - 88 ) 



Comparing (2-88) and (2-16) yields the scattering ampli- 
tide* 



*In order that the comparison be made the correct form 
of G(r,r') must be determined. 
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V(r ' ) dr ' 



(2-89) 



f(O-) = 



= 2mf iq.r 

' fi 2 i 



■* rV 



where q = k^ - k^. For a spherically symmetric potential 
V, (2-89) is 



f(O-) = 



2m 

<2 2 
ft q 



O 

1 



r ’ sin qr 1 V(r ' ) dr' 



( 2 - 90 ) 



For a point charge of Ze, the potential V is Ze /r, thus 



f pt (S ' ) = 



(2-91) 



while for a distributed charge of Ze of the form 
2 



V D " 



Ze 



(1- e" ar ), a ^ 0 , 



2 2 
^ _ 2mZe a 

.22 2 , 
ft q q + a 



2 - f pt. (W 



a 



2,2 

q + a 



( 2 - 92 ) 



Hence the structure of the distributed charge shows up as a 

factor times the point charge scattering amplitude. It can 

be shown that the form factor here is the Fourier transform 

of the charge distribution. Figure 6 illustrates the 

2 

behavior of (2-91) and (2-92) as a function of q . The 

p 

slope of f ( 0-) at q = 0 is proportional to the radius of 
the structure under consideration. Before continuing it 
must be stressed that all that has been said so far of form 
factors refers to the non-relativistic case. In a rela- 
tivistic example difficulties are immediately apparent. A 
charge distribution, spherical in one reference frame is 
elongated in another frame so that the meaning of the radius 
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Figure 6 

of a charge distribution is not clear. No longer can a form 
factor be thought of properly as the Fourier transform of 
the distribution. One has to be concerned with the form 
factors themselves and not their three-dimensional Fourier 
transforms in a particular reference frame. 

The Meson theory of nuclear forces mentioned briefly in 
Sections I part D and II part A3 states that virtual pions 
are constantly being emitted and reabsorbed by a nucleon. 
Further, two nucleons are stuck together through an exchange 
of pions. Now, since nucleons are surrounded by a cloud of 
pions, the virtual photon emitted by an electron scattered 
from the nucleon might have in fact, interacted with the 
nucleon meson cloud. So the simple electromagnetic inter- 
action of a photon and a nucleon may be quite a bit more 
complex and depend on the dynamics of the strong interaction 
of the pion and the nucleon. Dispersion relations then 
attempt to understand form factor behavior as a function of 
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2 

q in terms of the strong interaction dynamics. What 
follows is a simple minded attempt to make plausible the 
form of the equations put forth for the form factor. 

The basic idea of dispersion relations (spectral repre- 
sentations , or S-matrix theory of strong interactions) is 
that a quantum mechanical amplitude for a physical process 
is the boundary value of an analytic function of a complex 
variable. For electron scattering F(q ) would be the 
boundary value along the real line of the analytic function 
F(z) where z -* q + ie. A dispersion relation is an inte- 
gral equation that relates a dispersive process to an 
absorptive process. The best known example is that of the 
dispersion of light in a dielectric, where the complex index 
of refraction n(to) = n (to) + i rr- a(to) is expressed as an 
integral over all frequencies involving the absorption 
coefficient a(co) : 



n ( w) - 1 = — p/~ g - ’£ m 

r ir / • 2 2 

J n to - to 




a (to' )dto r 
.2 2 



to -to 



( 2 - 93 ) 



(P represents the Cauchy principal value of the integral) 



A classical optical example of ( 2 - 93 ) is Sellmiers equation 
for the anomalous dispersion of light. 



2 



n 



1 - I 




( 2 - 9*0 



where A is the incident light wave length, A^ the i th 
resonant wavelength, A^ a constant, and N the number of 
natural frequencies of the medium. 
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In electron scattering 



F(q 2 ) 



1 

IT 



00 

I. 



ImP(q ) 
Hm 2 'J' 2 -'! 2 




(2-95) 



TT 

is a direct analog of (2-93). Now three mesons, the p, <p , 
and the u-meson have been found experimentally to have the 
correct properties so that they may couple to the nucleon 
and to the virtual photon. They correspond to the N natural 
frequencies above. Hence it is plausible to write keeping 
(2-9*0 in mind 



'(q 2 ) - j 



j = P, n. 



( 2 - 96 ) 



i= 1 1-^2 

m. 

l 



In terms of the physical form factors G, (2-96) might be 
written 



G e (q j = 



®wp g 4>p ^pp 



and 



2 + 2 + 2 + 1 g wp " s *p ‘ s pp 

1 -h '-h '-h 

m rn 

0) <J) p 



G ( Q 2 ) = + . y <t>P + P .P... + i , u _u _u (2-97) 

M'‘ qj 2 2 2 y iop y <j>p y pp u 

x 2 2 2 

m m . in 

w <}> p 



where the g's and y's are coupling constants to be fitted to 
experimental data. The twelve coupling constants in (2-97) 
are not all independent and their number can be significantly 
reduced. The best fit available for the proton was used 
here to predict G^ for use in (2-8*1). De Vries [28] adjusted 
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the free parameters appearing in (2-97) so that a minimum 

2 

is obtained in a x fit to experimental data. The following 
values were determined: g^ = 2.628, g^ = 1.853, 

gpp = 1.191, where = 15.6 F - ^, m^ = 26.6 F~^ are 

experimentally confirmed resonances of fairly narrow width. 
Since the p resonance appears as a rather broad peak (750 ± 
100) MeV and since some observers [29] have shown that the 
effective position of the p mass can be at about 600 MeV, 
the mass of the p resonance has been taken as a free 
parameter, determined to be m^ = 8.^63 F . Before con- 
tinuing it should be noted that both (2-9^0 and (2-97) are 
approximations of (2-93) and (2-95) respectively in the 
sense that the resonances considered are taken to be delta 
functions in the integrands of the more accurate form factor 
expression (2-93) and (2-95) 

It is interesting to note that the theory outlined above 
actually preceded the discovery of the mesons that it 
requires. Dispersion relations as applied to the strong 
interaction was developed in 1957 by Chew [30] et al on the 
basis of a two-pion exchange approximation. This theory 
was not able to explain form factor behavior, and it later 
showed (1959) that a three pion resonance must exist in 
order that theory and experiment agree. 

A multiplicity of form factor notations exist today. 
Coupled with their often misleading names and differently 
defined static limits, they present a baffling array. 
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Commonly, the starting place is with the two functions 
F -^ ( q ) and F 2 (q ), called the Dirac and Pauli form factors 
(really there are four form factors, two for the proton, 
two for the neutron) . They describe the structure of the 
normal Dirac coupling and the anomalous Pauli coupling 
between the nucleon and the electromagnetic field. and 

F 2 are not directly related to the electric and magnetic 
moment distributions of the nucleon, but in older papers 
they are often erroneously called electric and magnetic 
form factors. Still another general name for F^ and F 2 
are vertex form factors referring obviously to their use 
in the current operator J^. The Dirac and Pauli form fac- 
tors have the following static limits: 

F 1 (0) =1 F 1 (0) = 0 

P n 

F 2 (0) = 1.79 = k F (0) = -1.91 = K n (2-98) 
p p n 



Hand, Miller and Wilson [27] popularized the use of the 
physical form factors G. They were introduced to simplify 
the Rosenbluth cross section by removing the cross product 
terms in F^ and F 2 and are in no way more fundamental than 
the vertex form factors. They are defined 



G E " F 1 + nlcF 2 _ P CH 



G m " F 1 + kF 2 ~ ^ 1+k ^ F MAG 



(2-99) 



+2 2 

where n = -q /^Im , m the neutron or proton mass, < the 

anomalous proton or neutron moment. F , and F m „^ were used 
^ ch mag 
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briefly before the G factors became popular. The form 



factors Gg and G rn do merit the names, electric and magnetic 
form factor, and they have static limits 



It is often found that still another type of form factor 

is convenient to work with. In isotopic spin notation 

where x = 1 for the proton and x = -1 for the neutron, the 
z z 

charge and magnetic form factors can be written succinctly 
for both the neutron and proton as 





G m (0) = 2.793 



G (0) = -1.91 
m 



(2-100) 



P 



n 




n 



n 



or as 




Now define the isoscalar form factors as 




1 




(2-101) 



and the isovector form factors as 



62 



(2-102) 




2 ^ °E 



G E ^ 



n 



1 

2 



G V = i(G 
m 2 m 



- ) 



m 



2.35 



n 



These (iso) form factors are convenient for use in theoreti- 
cal work because the electromagnetic current may be split 

s V 

into isoscalar and isovector parts, i.e. J = J + J 

MM M 

The intermediate states in the photon-nucleon vertex can 

be characterized by this isoscalar-isovector idea in that 

the two pion resonances with (T=0, J=l) to and <f> manifest 

s s 

themselves in the isoscalar form factors G^ and G , 

E m 

whereas the third resonance (T=l, J=l)p contributes to the 

isovector form factors G^ and G^. It is interesting to 

note that only the isoscalar form factor 2G„ = (G^ + G^ ) 

P n 

appears in the deuteron cross section, which might be expec- 
ted since the deuteron current is an isoscalar. 

There is an important relationship between the physical 

form factors that was discovered in 1961 at Orsay by 

2 

Lehmann, Taylor, and Wilson [31] • In low q measurements 

2 — 2 

(actually the work was at q = 2.98 F ) of the e-p cross 

section the following was noted: 

G m 

n frl 

G E = — 2- ( 2 - 103 ) 



_2 

Later data verified this up to approximately 100 F , and 
(2-103) came to be known and called the "scaling rule." 
There is some suggestion of (2-103) when the vertex form 
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factors are expressed in terms of the physical form factors, 
for example inverting (2-99) yields 



G — r)G G — G-r-, 

p = _E m p = _m E 

l l - n 5 2 l- n * 

The above have poles at n=l unless (2-103) is true. Drickey 
and Hand have experimentally verified (2-103) to include 
the deuteron with 



G„ m d 
E , - 



M 



d y' d 5 y d’ 1,71 “ M_ y d' 



(2-104) 



A further relationship of this type has been noted, 



G 



m 



G 



J2. = 



m 



n 



y 



n 



X 



(2-105) 



Figure 7 shows that the scaling rules can be used as an 

_2 

established empirical fact (to at least 20 F ) without 
introducing appreciable error. 
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Figure 7 



3 . Elastic Electron Deuteron Scattering 



Following the same general procedure as Rosenbluth 

did with the proton, Jankus (1956) derived [3] the cross 

section for elastic electron-deuteron (e-d) scattering from 

a point deuteron in the first Born approximation. His 

results differ qualitatively from the proton cross section 

in one important way, i.e., there are three form factors 

characterizing the deuteron while two suffice for the proton 

This is in agreement with the work of Glaser and Jaksic [32] 

Their study showed that the cross section for scattering a 

relativistic electron from a potential with spin J, contains 

2J + 1 form factors. The three deuteron form factors are 

found to correspond to the charge, quadrupole, and magnetic 

2 

moments in the static limit q -*■ 0 . The general form of the 
cross section is as it should be for the one photon exchange 

2 2 2 6- 

process, A(q ) + B(q ) tan ^ • 

Specifically Jankus found that 



da 

dft 



where 



^d^ /: 



n , 2E . 2 ©■ 
1 + ^ Sln 2 



i-l 



Mott \ M^c 



D^(q 2 ) + | n 2 D 2 (q 2 ) 



+ |n { 1 + 2U-n) tan 2 |) D 2 (q 2 ) 



( 2 - 106 ) 



V q) = / [ u 2 ( r ) + w 2 (r)] j o (x)dr — 1 (2-107) 

J o 

is the charge structure factor. 
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„ ,+2v 6/2 m 2 

d q q ’ * — M d 



°o ^ 

y [u(r)w(r) - j 2 (x) 



dr 



M^Q (2-108) 



is the quadrupole moment structure factor. 

CO oo 

dJJ(5 2 ) =J [u 2 (r) - |w 2 (r) j j 0 (x)dr + ~J |u(r)w(r) 



+ — j/ ~ | J ? (T)dr -1 - h 

✓ 2 J 2 



D 



(2-109a) 



accounts for the contribution of the intrinsic magnetic 
moments of the proton and the neutron to the scattering 
process and 

OO 

D^(q 2 ) = |- f w 2 (r)[j o (T) + 0 2 ( T )] dr * \ P p (2-109b) 

o 



is the magnetic contribution to the scattering process 
arising from the convection of charge in the deuteron. 
Together (2-109a,b) are 

d m<5 2 > = + V 2D m + d m] <2 - 110) 

p 

Also given are the static limits (q -»-0). The j and are 
spherical Bessel functions, t = |q|r/2, 0 is the quadrupole 
moment and the percentage D-state (see page 127). 

While the currents in e-p scattering were easily found, 
the deuteron current presents somewhat more of a problem. 

The virtual photon is expected to interact with one nucleon 
only. But the interaction amplitude of an electron with a 
bound nucleon is unknown. The natural thing to do is to 
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assume that since the deuteron is a loosely bound system, 
the photon interacts with a nearly free nucleon, while the 
other nucleon participates just as a spectator. This is 
called the impulse approximation and is described graphic- 
ally by Figure 8. 






Here the bubbles represent the photon - free nucleon vertex. 
Distortion of the free electromagnetic form factor (vertex) 
due to binding, i.e., the difference between bound and free 
form factors, is thought to be less than 5% [33]* Now that 
the impulse approximation has freed the nucleons, the nucleon 
participating with the photon in the interaction is described 
by the free field Dirac Equation. 

Since the impulse approximation involves the inter- 
action of a photon with an individual "free" nucleon, the 
current is again of the form UO^U, where O^ 1 is given by 
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( 2 - 111 ) 



°p - F lp«! 2 >V U + P 2p (q2) ° m ’'V 
°n " P ln (q2) ^ + P 2n< q2 > qWVq V 



The deuteron current is then 



r P 



-uo^u + vojjv. 



(2-112) 



Adler [45) writes the relevant part of the S-matrix form in 
the impulse approximation as 



S fi | 2 - |M fl | 2 - Up ^l 2 

2 . 



(2-113) 



Just as in the proton case, is evaluated and found 

to be 

l M fi! 2 = G c + jp 2 °Q + | nG M [ 1+2 ( 1 -n) tan 2 f](2-ll4) 

where p = A -q /4M d , and G c , Gq and G^ are the charge, quad- 
rupole and magnetic form factors given by (including their 
static limits ) . 



G c 1 ^ G E + G E ^ D c 2 n 
p n q = 0 



(2-115) 



and 



g q ^ g e + g e ^ V 

p n 



M a Q 



(2-116) 



and 



g m ^ g e + g e + ^ g m + g m > 2d m 

p n p n 



^ , 

M y d’ 
P 



(2-117) 
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where the D's are given in equations (107) through (110). 

The steps from a point deuteron to a finite size deuteron 
is thus carried out by the introduction of the free nucleon 
form factors. This step is justified by the impulse 
approximation . 

2 2 

The contributions of (charge), (quadrupole) and 

2 

G M (magnetic) to 

A(Q 2 ) = G 2 (q 2 ) + | n 2 G 2 (q 2 ) + |nG 2 (q 2 ) (2-118) 

are shown in Figure 9. 

k . Neutron Charge Form Factor from e-d Scattering 

Equation (2-117) can be simplified considerably by 
scaling, i . e . , 

g m h = y d<V ’ y d = 

d d 

and 



°M °M - 
— E. = _a = °E 

"p v n P 



For the case of the deuteron 
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(2-119) 



(2-120) 
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Figure 9 

2 8 2 2 2 2 
The functions (charge), — n Gq (quadrupole) and nG^ 

(magnetic). The smooth curve is the sum of these three 

terms (taken from Ref. 46). 
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2 

where is called the deuteron structure factor and is the 
same as the result of Jankus for a point (nucleons) in the 
deuteron. From (2-120) 




The aim of the experiment is the evaluation of the 

charge form factor of the neutron. The magnetic contribu- 
2 2 r 

tion to G. and G must therefore be eliminated. In the case 
d p 

of the deuteron the magnetic contribution is given by 
(2-120) as 



n 2 1 2 / , , « . 2 G-v 

C = - nu d (1+2 tan -) 



so that 



G 
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E 



d 



G EXP 
1 + C 



with 




c 

exp 

a mott 



d 



(2-122) 



(2-123) 



where the quadrupole term has been neglected. 

In the case of the proton the magnetic contribution is 
found by considering 



G 2 . G m 

0 E O' M 0 

G 2 = —2. + n t- 2- (1+2 tan 2 f) (2-12^) 

P l-n 1-T1 2 

Scaling assumes G^ = yG E where y = 2.79 so, 

°E 

G 2 = |l - ny 2 (1+2 tan 2 — )| . (2-125) 



The term 

C = ny 2 (1+2 tan 2 



( 2 - 126 ) 
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is then the magnetic correction term. Thus 



G = 



E 

T* (l + 

l-n 



c) 



so that 



G 




G 2 (l-n) 

p y 

1 + C 



(2-127) 



where 

q 2 _ exp 

P a mott 

Finally, for experimental comparisons 




or 



°E . G E. 
n _ 1 d _ , 

G E " P d G E " 



( 2 - 128 ) 



The value of G^ from (2-128) is very sensitive to the 

n 

uncertainty in F^ and G E /G E . For example a 1 % error in 

n p 

F^ translates into about a 10# error in G^ /Gg . 

n p 

The only reliable data in the range of considered 
here have been presented by Drickey and Hand [1] . Their 
results were in complete disagreement with the neutron- 
electron interaction slope. Their data will be presented 
after reanalysis of the theory in part B5 . 
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5 • Review of Relativistic Deuteron Theories 



Essentially , all the material above is based upon 

the solutions to a nonrelati vistic Schroedinger equation for 

the two-nucleon problem. The results produced using the 

nonrelativist ic S and D state wave functions are, at times, 

surprisingly good. Still, there are problems as the nagging 

difficulty in the determination of the percentage D state, 

and the theoretical determination of the deuteron magnetic 

moment, not to mention the conflict that is the point of 

this work, i.e., G p * 0 while dG-. /d(q 2 ) 

n E n h - 0 • 

Hence sufficient reason exists for improving the theory, 

as well as the experiment . 

The Jankus theory on which the cross section (2-106) 
was derived is itself one of the major limitations in the 
development of better answers for the deuteron. Gross 
states that the theory gives no indication of what has been 
left out, or how to calculate corrections, which might 
include the effects of of f— the— mass— shell contributions 
arising from the fact that the nucleons are bound and not 
free. Intermediate states, i.e., pion exchange currents, 
could cause further corrections. Thus to isolate possible 
corrections a fully relativistic theory should be developed, 
one that would reduce to Jankus theory at low q 2 . Correc- 
tions then could be localized and more carefully defined. 

This is a difficult order to fill, and to date no entirely 
satisfactory solution has been proposed, though many authors 
have tried. The most interesting attempts are briefly dis- 
cussed below. 
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Van [3^] has calculated, using deuteron wave functions 
derived from a Bethe-Salpeter equation, relativistic cor- 
rections to the deuteron form factor. These wave functions, 
which are relativistic generalizations of Hulthen wave 
functions (no hard core) tend to make the functional D of 
Equation (2-106) larger at high momentum transfers than do 
models having a hard core. The essential part of Van's 
results is a comparison between his predictions of (F^ 
in Van's notation) calculated with the (relativistic) 

Hulthen wave functions and the corresponding D ( -,(F^, Van) 
of the nonrelativistic case. The comparison showed 
^ch (3^ > -^ch a ran £ e 9^ from 0 to 9 F - ^*, which 
results in F^ being larger than in the relativistic case. 
Recalling (2-128) quickly confirms the fact that the end 

result is to make Gg /Gy even smaller. The increase in 
ch NR n P 

F c ^ from F ch is attributed to the quadrupole contribution 
2 -2 

(+3.5$ at q = 1 F ), and to final state wave function 

2 -2 

Lorentz contraction (-1.5 % at q = 1 F ). 

Adler [^51 has obtained what amounts to first order 
relativistic corrections to the nonrelativistic impulse 
approximation deuteron form factors. His theory is intrin- 
sically nonrelativistic and differs from that of Section 
II part B2 only by the use of Pauli spinors (for convenience) 
instead of Dirac spinors , and by retaining powers of 

deuteron moments such that the final cross sections are 

2 2 

correct to order q /m , The extra terms now found in the 



7 ^ 



deuteron current are incorporated into additive corrections 



to and Gq. Specifically 
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Q 



* 2 \f 
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5 / 2 - 4 , 



->2 2 
q r 



dr 



q|r 



and where 



g rs = G r + G r , G R = P 1 (q 2 ) + 2KF 2 (q 2 ). 
p n 

At q 2 = 1.0 F -2 where G c = 0.55, AG^ - 8 x 10 -5 and where 
R _2| 

Gq = 0.03, Gq = 2 x 10 , hence the corrections are rather 

too small to effect F^ greatly but go in the right direction. 
Further, these corrections to F, do not alter the value of 
y^ since it depends only on P R , and the value of is 
changed by a bare 1% . In summary then, it is doubtful that 
these corrections are very meaningful in a theory using 
nonrelati vistic wave functions. 

Adler glossed over at least one sticky problem in the 
above work. He admits the need for a "fudge factor" in 
resolving the kinematical difficulty in the lab frame of an 
energy transfer (recoil) to only one of the nucleons by 
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o 



the virtual photon. Gourdin [35] has made what amounts to 
the same calculation but in the Breit frame where the 



nonrelativistic deuteron wave functions into the Breit 
frame. His solution to this was to assume that the wave 
function were unchanged. Gourdin' s results were approxi- 
mately the inverse of Adler's, i.e., of the same magnitude 
but opposite in sign. Unlike Adler, Gourdin did derive a 

very small correction to the magnetic form factor but it is 

2 -2 

so small that it is not observed for q £ 6 P . The value 
then of these works lies in the fact that it almost certainly 
seals the fate of any straightforward attempt to determine 
(quasi) relativistic corrections. 

The strictly relativistic approaches (a covariant 
dispersion relation representation) are more basic and more 
difficult than the quasi-relativist ic work above. Success 
with this technique is limited to the low momentum transfer 
region where it is now believed that the theory is more 
accurate than experiment. The early papers of Jones [*l] and 
Cutkosky [36] are interesting when applied to the Gp 

n 

problem. Using the triangle impulse approximation diagram 
(page 67) Cutkosky calculated the simplest absorptive 
amplitude which yields a deuteron structure factor of (all 
spins ignored) 



2 -*2 

momentum transfer is entirely spacelike, i.e., q = -q , 
2 

q^ = 0. His difficulty was in the transformation of the 




tan — ,a =/wM. (2-129) 




/l+n 
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The relativistic analog using the tail approximation is 



F^(q 2 ) = / u ( r ) j q ( qr/2 ) dr 

J o 



where 




tan 




( 2 - 130 ) 



a 



a Mott 





(2-131) 



T A 

Comparison of the two F^ and F^ yield the relativistic 
corrections 



and 





l+n 



/1+r) 



( 2 - 132 ) 

(2-133) 



These imply that a surprisingly rapid curvature is needed at 
2 

small q for Gg in order that the n-e interaction be in 

n -2 
consonance with Gg - 0 at 0.3 - 0.6 P 

n 

Gross has done the most extensive work in the area of 
the purely relativistic approach [2 , 8, 9, 10, 11]. He 
ambitiously set out to construct a complete theory fully 
relativistic and fundamental in the sense that very few 
phenomenological parameters occur and based on single- 
variable unsubtracted dispersion relations and coupled 
unitarity equations. He had partial success but numerical 

predictions of the ratio of relativistic to nonrelativistic 

2 

form factors seems much too large at high q . In the range 
2 —2 

of low q (=0.25-1. OF ) agreement is reasonable, and 
his predictions for y^ and are within 2 % of the 
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experimental value. In work on relativistic corrections to 
the impulse approximation in elastic e-d scattering Gross 
realized that the purely dispersion technique tended to 
mix up the internal structure of the interacting particles, 
i.e., the deuteron wave function and the nucleon form fac- 
tors with the deuteron current. This leads him to eventually 
express his theory in a formulation closer to that of 
potential theory. 

Specifically he obtained a modified Bethe-Salpeter type 

-2 

wave function retaining terms up to the order of M , and 
then used the relativistic nucleon current to derive the 
deuteron form factor. An estimate of the size of the cor- 
rection to the charge form factor is obtained by taking the 

p 

slope of the correction term at q =0. Gross found that 

2 

the correction is proportional to -1/8M which is of 
right sign so as to make the neutron charge form factor 
more positive (towards agreement with thermal n-e slope 
data) . 

The relativistic correction mentioned above came mainly 
from relativistic modifications of the deuteron wave func- 
tions and the nucleon current. Adler had previously looked 
into the nucleon current modifications but his analysis had 
been based on nonrelativistic deuteron wave functions. 

Gross states that the wave function modifications have the 
largest effect and hence Adler’s corrections are insufficient. 
Van did use relativistic wave functions from a Bethe-Salpeter 
equation in his study of corrections to the impulse 
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approximation. His difficulty was in using a separable 
potential which was equivalent to using a Hulthen wave 
function which has no hard core, and thus gives values of 
too large. The trouble is that here relativistic cor- 
rections to e-d scattering and the model dependence 
question have been intermixed. Van did not avoid this by 
comparison of his relativistic corrections to corresponding 
nonrelativistic models since the two were not in the same 
reference frame. That is to say the relativistic model (in 
the Breit frame for example) when Lorentz transformed into 
the deuteron rest frame was not compared to a corresponding 
nonrelativistic model, which if it is a good model already, 
simulates all relativistic effects as seen in the rest 
frame . 

The principal contribution to the correction of the 
charge form factor of the deuteron comes from taking into 
account the recoil motion of the deuteron, which in turn 
distorts the wave function. Lorentz contraction also adds 
to the wave function distortion, though not as much as in 
Van's case. Ignoring the correction for the quadrupole and 
magnetic terms, the result of Gross' calculation for the 
charge form factor is 

G C = (G E + g E ) D C (2-13*0 

P n 



Van did this in order to improve on Jones and Cutkosky's 

e -ar 

asymptotic "tail" wave function approximation = . 
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where 



D R = (1+|) J (u 2 + w 2 ) j Q dr + AD C (2-135) 

0 



and 



A°c = 



oc 

• 1 / 



2 2 

(u - uru' - uku + 7w - wrw ' - wkw) ( J q + ^ dr 



oc 

tjf 



(uru ' + wrw' ) j dr, 



(2-136) 



In this expression prime refers to differentiation with 
respect to r, and 



k = r 



2 





(2-137) 



The arguments of the Bessel functions are of course qr/2, 
and q is | q^ | . Now in terms of Gross's corrections, A(q 2 ) is 



A(Q 2 ) = (G +G e ) 2 d£ 2 = (G e +G e ) 2 (F 2 +AP d ) (2-138) 
p n p n 



where the quadrupole contribution to the charge form factor, 
8/9n Gq, is small. Thus, since A is an experimental quan- 
tity and does not change due to corrections, a change is 

forced in Gg . To see this write 
n 

G c = (G e + Gf)Df = (G e + G e R )D c (2-139) 

p n p J n 

(no corrections). 
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But 



G c = CG e + G R )D R = (G + [G™ + AG ])(D RR + AD ) (2-140) 
p n p n n 

(with relativistic corrections). 

Equating the right hand sides of (2-139) and (2-140) yields 

-AGg D c - (G e + G e ) AD C - G e AD c 
n p n p 



or 



AG 



E 



n 



AD, 



-G 



E r> n NR 
P D c 




AF 



d 



F 



d 



(2-141) 



where F^ is the deuteron's charge structure function. Gross 

wrote (2-l4l) leaving out the F^. Since F^ 1 for very low 
2 

q , AG e - Gg AF^ gives an upper limit to the correction as 
2 n p 

q 4 0 in the n-e interaction. Figure 10 summarizes the 

2 

corrections. In the range of q considered here. Gross' 

2 

correction AG E is proportional to - and practically 

a 8M ^ , .. ...i-w 6 Cf v ^ ' 

model independent. p 

Using this correction to the charge form factor of the 

neutron plus Feshbach-Lomon wave functions for the deuteron. 

Gross reanalyzed the G E - 0 data of Drickey and Hand. 

n 

Figure 11 shows a comparison of results with Partovi (i.e. 
HJ) and Feshbach-Lomon wave functions. Gross' correction 
coupled with FL wave functions seems to remove the discrep- 
ancy between the electron scattering and n-e interaction 
results . 
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Figure 10 



REL 



The relativistic correction is labeled by AG„ , the 

b n 



difference between the Feshbach-Lomon mpdel and the Partovi 
model is labeled Gl^L). g£ p ) and the meson exchange contri- 



^n x 1Jn 

butions are labeled AG F 

^n 



The thermal neutron slope is 
given for reference (taken from Ref. 2). 
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(a) 



(b) 



Cc) 




The data of Drickey and Hand. The dashed lines are the 
thermal neutron slope, +0.0193 ± 0.0004 and the fit of 
Chilton and Urhane (unpublished). The number on the right 
of each figure is the slope determined from the first three 
data points of each graph. (a) The data analyzed with 
Partovi wave function, (b) the data analyzed with Partovi 
wave function including relativistic corrections, and 
(c) the data analyzed with the Feshbach-Lomon wave function 
including relativistic corrections (taken from Ref. 2). 
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A recent attempt by Chilton and Urhane [37] to explain 

G„ ^ 0 for q 2 - 0.3 - 1.0 F -2 with Gg (0) / 0, was made 
n n 

by fitting the nucleon form factors with rational functions 

2 2 

of the variables 1 + q /m^ . They found Gg is never far 

n 2 

from zero and bends over rapidly with increasing q (see 

Figure 10). Nevertheless with Chilton and Urhane’s fit, 

_2 

G p does not get back to zero till beyond 12 F 
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III. APPARATUS AND EXPERIMENTAL PROCEDURES 



A. THE NPS LINAC 

The Naval Postgraduate School linear accelerator (LINAC) 
is a low Intensity (20 yamp maximum average current), 100 MeV, 
three section, 30 foot long electron linear accelerator. 
While its construction and general operation have been 
described elsewhere [38] , several important modifications 
have been accomplished and will be mentioned briefly. 

A supplementary cooling water system with a steam heat 
exchange was added for the accelerating sections. This was 
done so that the water temperature could be precisely con- 
trolled (±1°C) throughout a run lasting for an extended 
length of time. In the present experiment, runs of up to 
38 hours were made, and since the machine's energy stability 
is strongly dependent on the accelerator's changes in 
temperature, this is a critical point. 

Until recently the target chamber vacuum had been 
isolated from the accelerator sections by a three mil 
aluminum window 32 inches upstream from the target. This 
window caused considerable beam spread and by the time the 
target was illuminated by electrons, the beam spot was 
approximately 3/16 - 1/ M inch in diameter. Now, with the 
addition of a fore-diffusion vacuum pump combination in the 
target area, the entire accelerator, from electron gun 
(source) to target, is open. In fact with the spectrometer 
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coupled to the target chamber, electrons pass only through 
the target Itself and a five mil aluminum exit window from 
the spectrometer before striking the detectors. Upstream 
secondary emission monitors (SEM) have been moved to loca- 
tions downstream of the target. The accelerator - target 
chamber - spectrometer combination is maintained at a 
pressure of 10 mm of Hg. The removal of the SEM and 
window has significantly reduced the beam spot size, and 
it is now estimated from beam appearance on a zinc-sulfide 
(ZnS) screen that 2/3 of the electrons are confined to a 
spot 1.0 mm in diameter. 

Past nuclear structure experiments [39] have experienced 
severe problems in background radiation levels. Major 
shielding modifications were made to correct this situation. 
More paraffin and borated paraffin in the form of a floor 
to ceiling wall was added to the existing shielding in the 
beam deflection area. In the target area, a platform loaded 
with paraffin four feet above the floor and extending 30 
inches out from the wall was suspended from the overhead. 

Its purpose was to shield the counter house from line-of- 
sight radiation from the beam dump. Wax was placed in the 
spaces between the I-beams in the overhead in order to 
reduce reflection of neutrons into the counting system. 

The primary shielding effort was on the counting house 
itself. Concrete blocks were replaced by boxes of paraffin 
and much lead was replaced by special bricks composed of 
5 % boron, 80 % lead by weight and polyethylene. A cross 
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section of the counter house is shown in Figure 12. This 
shielding arrangement reduced the background to quite an 
acceptable ratio of 2000 to 1. 

B. THE TEN-CHANNEL COUNTING SYSTEM 

1 , The Ten-Channel System 

In general there are two main categories of detec- 
tors: slit-defined momentum acceptance, and counter- 

geometry defined acceptance detectors. A single counter, 
slit defined detector was replaced by a ten counter 
geometry defined system. The ten counters, which span a 
portion of the focal plane of the spectrometer and define 
adjacent momentum intervals, are known as a "ladder." The 
advantages of this arrangement over the single channel 
system is obvious; for one spectrometer momentum setting, 
ten pieces of information instead of one are gathered. 

The ladder is composed to ten 7/16 inch high, 1/16 inch 
thick NE102 plastic scintillators placed approximately 1/2 
inch behind the focal plane of the spectrometer as shown 
in Figure 12. A large one-piece scintillator used in 
coincidence with the front counters and called the backing 
counter, is placed 3 inches behind the 10 singles counters. 
The choice of 7/16 inch high (in the dispersion direction) 
counters produce a momentum acceptance Ap/p of about 0.3$ 

*The ten scintillators (counters or channels) are num- 
bered in ascending order 1 to 10, where channel 1 represents 
the lowest momentum accepted. 
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per channel or about 3 % for the entire ladder. An electron 
scattered from the target and subject to the spectrometer’s 
magnetic field B, follows a path of radius r, given by 




Thus higher momentum electrons are seen by the higher 
numbered channels. Further the momentum resolution is 
greatest for the higher channels. An electron traversing 
the spectrometer and passing through a particular scintil- 
lator produces a pulse that is collected by an Amperex 
XP 1110 photomultiplier tube, amplified, transmitted, 
discriminated, fanned out, and finally registered on one 
of ten scalers. A block diagram of the counting system is 
shown in Figure 13.* With the exception of the scalers, 
which allow a counting speed of 20 MHz, the electronics is 
capable of 100 MHz counting rates. 

2 . Resolution and Efficiency of the Ten-Channel System 
In order that the system be useful, two quantities 
must be known: The relative counting efficiency of each 

channel and the momentum separation of adjacent channels. 

The momentum separation (resolution) of the channels was 
obtained by taking elastic scattering peaks from a thin 
carbon foil and detecting the electrons scattered as a 

*The Chronetics and EG and G counting equipment was 
loaned to NPS by the High Energy Physics Laboratory, Stanford 
University . 
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function of spectrometer momentum in each channel. That 
is, the 10 counters were treated as if they were single 
counters while the sharp elastic peak was moved through 
them by changing the momentum setting of the spectrometer. 

A plot of counts in the scalers vs. spectrometer momentum 
yields 10 elastic scattering peaks (Figure 14). The separ- 
ation of the peaks as a function of electron momentum p is 
a measure of the channel's resolution. The 10 elastic peaks 
produced are not identical in height or area, but since 
they all "see" the same elastic carbon peak, the differences 
are ascribed to the differing relative efficiencies in detec- 
tion and recording of the electrons. 

Now specifically, the resolution of the 10 channels was 

determmined by taking 21 - ten peak sets at incident elec- 

MeV i MeV 

tron momentums from 27 to 9^ and with scattering 

c c 

angles varying between ^5° and 105°. These experimental 
data were fit to a quadratic expression 

p = p ( a + BN + YN 2 ) (3-2) 

o 

where N is 'the channel number, p is the momentum setting of 
the spectrometer, and p is the mean momentum of channel N. 

A best fit produced values for a, B and Y. If p is defined 

o 

as the momentum setting for which a particle of momentum p 

s 

entering the spectrometer passes through the center of 

channel 6, then (3-2) can be rewritten as 

■ 7- \ = P [1 + A(N-6 ) + B(N-6 ) 2 ] 

p N (p s } S 
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scattered electron energy 
(MeV) 

Figure 14 



A and B again are constants independent of p. The channel 
momentum acceptance is then 



d( PN (p s >> / 

d(N-6 ) / P s 



A + 2B(N-6) . 



(3-4) 



Typical values of A and B are 3.0 x 10 - ^ and -1.3 x 10 - ^ 
respectively and it is clearly seen that the higher channels 
(i.e., N > 5) have a higher resolution than the lower 



channels . 

So that the energy spectrum of the scattered electrons 
can be determined properly the relative efficiencies need 
to be measured. The efficiencies can vary for several 
reasons, some of which are: Differences in detection 

thresholds and physical properties of the plastic scintil- 
lators; spectrometer transmission characteristics; and 
solid angle. The determination of the relative efficiencies 
was based on the method of Crannel and Suelzle [40] , and 
essentially involved measuring the shape of an inelastic 
continuum of scattered electrons. The shape was found by 
observing the spectrum with all ten channels at several 
different momentum settings assuming initially that each 
channel has efficiency one. Then the counts from individual 
channels are compared to the calculated shape and corrected 
channel efficiencies are obtained. Next, a more accurate 
shape is predicted on the basis of the corrected efficiencies 
and the process is repeated. Three or four iterations will 
normally produce constant efficiencies. 
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Typically, five different momentum settings over a 
total momentum range of 6 % (twice the 10 channel span) with 
at least 1000 counts per channel per setting were taken. 

It is important that the efficiency data be taken in the 
energy range in which the experiment is to be performed. 
Hence if an elastic deuteron peak is located at 50 MeV, then 
the machine energy should be adjusted so that a smooth part 
of some inelastic continuum lies at that point. In carbon 
12, the radiation tail of the elastic peak is sufficiently 
smooth at 30 MeV below the peak, thus for the above situa- 
tion a machine energy of 80 MeV would be set, and the five 
momentum (energy) settings could be taken at 52, 51, 50, 

49, and 48 MeV for example. Carbon is obviously not suited 
for the NPS machine because then a practical limit of 
65 MeV would be imposed on the incident energy for the 
proton and deuteron peaks (i.e., 65 + 30 MeV = 95 MeV 
needed for efficiency and any more would exceed machine 
capability). Aluminum was substituted; its inelastic tail 
contains peaks only at 2.2 and 3.0 MeV. Arbitrarily a value 
of -6 MeV down from the aluminum-elastic peak was chosen as 
the experimental area from which to obtain the relative 
counting efficiencies. Normally it is best to exactly dup- 
licate in the efficiency measurement the conditions of the 
particular experiment, but this was found to be impossible 
in some cases. On the high energy runs, 94 MeV (location of 
elastic peaks under consideration) , 94 + 6 = 100 MeV was an 
impractical energy to obtain and maintain. Thus the 
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efficiency data were taken at a somewhat lower energy. Also 
aluminum not or was always used for efficiencies 

since it was desirable not to take the risk of burning or 
destroying the experimental targets during the long high 
beam strength used in efficiency runs. 

To calculate the relative efficiency of channel k it 
is necessary to predict the total number of counts channel 
k receives in the L spectrometer momentum settings. Thus 



p h (k .p s ) 



SW | / E k (p.p s )T(p s )dp (3-5) 

s=l s=l p £ (k,p s ) 

where p is the particle momentum, T(p )dp is the number of 

s 

electrons with momentum between p and p + dp which enter 
the spectrometer, E^(p,p g ) is the efficiency of channel k 
for particles of momentum p, and p h (k,p g ) and p^(k,p s ) 
the highest and lowest momentum accepted by channel k. 
Assume that 



T( p s ) = 1 V 



n+1 

s 



(3-6) 



n=0 



where R can be determined (needs to be done only once) by 

n 10 T e 

a least squares fit to the points T (p ) = £ C, (p ) as 

c s 1< = 1 ^ ^ 

a function of p . Hence the efficiency of channel k is 
given by the ratio 
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(3-7) 
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which is the total number of counts observed vs. the pre- 
dicted total number of counts. 

C . THE TARGETS 

Solid polyethylene (C^) an< ^ deuterized polyethylene 
(CD ) p were chosen as the sources of protons and deuterons. 
The immediate difficulty is that these plastic targets would 
be destroyed by heat (melted) in the concentrated electron 
beam. Hence the CH^ and CD^ targets were made to move in 
the horizontal direction so that the beam sampled a width 
of about one inch over the target face. The carbon in CH 2 
and CDg has much greater atomic weight (12 umu) than does 
either isotope of hydrogen (H = 1.00783 umu and D= 2.01^10 
umu) and hence suffers much less recoil when struck by an 
electron. Thus both the elastic proton and deuteron peaks 
ride well down on the inelastic tail of their respective 
carbon peaks. In order to isolate the proton and deuteron 
peaks it is necessary to subtract out the carbon content. 

A carbon target is needed then for standardization. Hence 
arrayed down the target ladder were movable CD^ and 
targets, a stationary carbon standardization target, alum- 
inum target for efficiencies and a ZnS screen for machine 
tuning and beam observation. 

The CD^ target was 3 cm x 10 cm x 0.10122 cm and 
specified to be 99 . 9 % enriched. The thickness of the CH^ 
and C targets were governed by the desire to make the number 
of C atoms in each of the three targets as nearly the same 
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as possible. With the material available a CU^ target 
0.10^3^ ± 0.00070 cm thick and a 0.09398 ± 0.00061 cm thick 
carbon target were constructed. The carbon target was 
produced by laminating 7 thin (5 mils) graphite foils. The 
densities of the samples were determined to be: p(CD 2 ) = 

1.081 ± .006 gm cm’^, pCCH^) = 0.925 ± 0.005 gm cm - ^, and 
p(C) = 1.028 ± 0.005 gm cm"^. 

In the course of a run the electron beam illuminated 
the CH^ and CD^ targets for as long as 8 hours apiece. 

After these long runs considerable browning of the material 
was evident, and the possibility of damage to the target 
because of a change in the composition of the material was 
considered. Subsequent runs on the undarkened area versus 
the darkened area gave no indication of any decomposition. 
Further the darkened areas tended to fade in time so no 
permanent damage could be supposed. The browning of the 
target is assumed to be an atomic electron dislocation 
phenomena . 

D. .COLLECTION OF DATA 

After deciding upon the momentum transfer to investigate, 
the appropriate incident energy-scattering angle combination 
was chosen. The basic measurement desired was the ratio of 
the electron-proton to the electron-deutron scattering 
cross section. Forward scattering angles were chosen where 
possible in order to emphasize the electric interaction. 

With this angular restriction the momentum transfers ranged 
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p _p 

for the proton and deuteron from about 0 . 1 <_ q < 0.5 P 
for the available machine energies. 

The machine was tuned for maximum beam and current 
stability. The energy defining slits were set for 
AE/E - 0.5%, and the beam spot was carefully steered to the 
target center line. Collection of raw data consists of 
recording the number of counts shown in scalers 1 through 
10, the backing scaler counts, the time elapsed while 
counting, the spectrometer setting in MeV, and the number 
of electrons passing through the target. All the scaler 
data are fed directly to a teletype machine that both punches 
an 8 track paper tape, and prints out a hard copy. The 
paper tape is then interpreted and IBM cards produced for 
use in computer data analysis. 

For this experiment data were taken in the following 
way. The machine was tuned to an energy approximately 6 MeV 
higher than the energy at which the elastic proton peak 
would occur. Efficiency data using a 20 mil aluminum target 
were taken. The machine energy was then lowered to the 
value selected for the proton run. Comparison elastic 
carbon peaks in CHg and C were taken to provide normaliza- 
tion data used in scaling the carbon for subtraction. Then 
the elastic proton peak from CE^ and the underlying carbon 
tail for subtraction were obtained. Machine energy was 
again lowered to the proper value and the deuteron runs were 
completed in a manner similar to that for the proton. On 
all peaks the counting rates in the front channels were 
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limited to 75 counts per second, while the backing channel 
was limited to 100 counts per sec. Tests showed that the 
counting system (presumably the scalers themselves, not 
the counting electronics) could digest counting rates of up 
to 200 counts per second before any appreciable loss was 
detected. Since the deuteron is twice as heavy as the 
proton, less incident energy is required to achieve the same 
momentum transfer as that of the proton. In all experiments 
a peak was determined by three spectrometer settings, hence 
30 points define a scattering peak. The small energy steps 
taken depended on the machine energy and the relative energy 
(momentum) separation of the various channels. The aim was 
to achieve a uniform distribution of points over the area 
investigated . 

Where the peak is located on the ladder is important. 

The end channels, 1, 2, 9, and 10 tend to have more widely 
varying efficiencies than the center channels and they 
experience edge effects, so they are unsuitable for the 
peak determination. It was found that placing the peak 
initially between channels 6 and 7 gave the best results. 
Then the top of the peak was defined by a single counter 
reducing a potential problem of matching the counts of 
different channels at a critical point. The peak height 
was chosen on the basis of time versus statistical accuracy. 
That is, a total of ^ x 10^ counts are needed under a peak 
to insure 0.5% counting statistics, but the time needed to 
accumulate that many counts is often prohibitive. 
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p 

particularly at the higher values of q . It was decided 
that no more than three hours per spectrometer setting 
would be allotted for counting. 

In order to extend the range of momentum transfers, 
another set of data was taken by the author at the Mark III 
accelerator at Stanford University. Where possible, over- 
lapping sets of data were taken to check for any systematic 

errors in the equipment. The range covered at Stanford 
2 -2 

was 0.20 £ q £ 0.80 F . No systematic deviations between 
the two sets of measurements were found and the data 
reduction was performed independent of the source. 

The experimental equipment at Stanford is basically the 
same as that at the Naval Postgraduate School with the 
exception that a 100 channel ladder is used and that part 
of the data reduction is performed on an IBM 7700 on-line 
computer . 
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IV. DATA REDUCTION 



A. ANALYSIS OF RAW DATA 

Raw data are accumulated In the form of counts in ten 
channels for a single spectrometer setting, and are analyzed 
or "unfolded" by matching the counts in a particular channel 
to the appropriate energy. The assignment of energy values 
to channels is determined by (3-3). The raw counts are 
first corrected for background, then counting rate correc- 
tions are applied. Finally the efficiency determined by 
(3-7) is used to predict a corrected count for the partic- 
ular channel. In practice all the above calculations are 
carried out by computer. Typical unfolded energy spectra 
for both proton and deuteron runs are shown in Figures 15, 

16 , 17, and 18. 

B. CALCULATION OF EXPERIMENTAL CROSS SECTIONS AND FORM 

FACTORS 

The experimental cross sections for e-p and e-d scatter- 
ing were determined from 



p,d 

a Exp 



N 



sc 



N. rn 

xn t 



2 -1 

-r-n K K, cm ster 
Aft s b 



(4-1) 



where N is the number of scattered electrons detected by 
sc 

the spectrometer, N is the number of incident electrons, 

in 

2 

n^ is the number of target atoms per cm , Aft is the solid 
angle subtended by the spectrometer entrance port, and 
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(MeV) 



finally K and K, are the radiative corrections of Tsai I4l] 

S O 

and Bethe-Ashkin [42], 

At the heart of the calculation of (^-1) is N which is 

sc 

the integrated area under the scattering peak in counts-MeV 
per MeV. The scattering peak is produced by drawing a smooth 
curve through the data points in the counts vs. energy plot. 

Both the elastic proton and deuteron peaks ride on the 
radiation tail of the elastic carbon peak (see Figure 15, 
etc.). Thus, in order to obtain the relevant areas for 
the cross section, the underlying carbon tail must be sub- 
tracted out. Comparison of the elastic carbon peaks in 
CH 2 and CD 2 to that of the pure carbon target provided the 
target normalization factors. These factors consisted of 
(1) "shifts" of the elastic carbon peak from the carbon 
target such that its peak energy coincided with that of the 
carbon peak from CH 2 or CD 2 ; and (2) vertical "scaling" of 
the carbon tail through a comparison of the areas of the 
elastic carbon peaks. The area remaining after the sub- 
traction is N 

sc 

The number of electrons incident upon the target is 
measured by a Secondary Emission Monitor (SEM) downstream 
from the target. The SEM is not an absolute device to 
measure the charge. In a ratio experiment, the only 
requirement is that it be stable during a set of runs. 

This stability has been frequently checked against a 
Faraday cup. Obviously the number measured is wrong by at 
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least the number of electrons scattered by the target, but 



this Is so small compared to as to be totally Insig- 

nificant. n^ Is a static parameter; It depends only on 
target composition and target angle <f>, 
pt N 



_ Q 

where p is the target density in gm cm , t thickness in cm, 
N Avogadro's number, A the atomic weight, and <f> the angle 

Cl 

between the beam line and the downstream target normal. The 

solid angle AQ is defined by 1 inch thick tungsten slits on 

the spectrometer entrance port, 16 inches from target center. 

It can be seen by comparison of two experimental cross 

sections that the only remaining quantities that enter a 

ratio experiment are N and n, . 

set 

As mentioned above, we have used an expression due to 
Tsai to account for the radiative correction K g . This cor- 
rection is a significant improvement over that given by 
Schwinger [ it 3 ] - Presented here is Schwinger's much simpler 
expression (*t-3) as an illustration of the parameters 
entering the correction. 
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Briefly, K is a multiplicative correction to the cross 

o 

section which accounts for radiation emitted during a large- 
angle scattering event (nuclear bremsstrahlung) . corrects 

for the radiation that occurs before or after the large-angle 
event and is sometimes known as thick target bremsstrahlung. 
For a discussion of the terms in (4-3) and (4-4) see Gordon 
[44] . The £ used above corresponds to the nucleon recoil 
factor in ( 2-77 ) . 

The ratio of the experimental cross section to the cross 
section predicted for a point nucleon (Mott cross section) 

p 

Equation (2-76), yields the experimental form factors G 

2 P 
or G^. To these quantities, the magnetic corrections as 

outlined by equations (2-122) and (2-126) are then applied. 

This yields the square of the electric form factor of 
2 

the proton, Gg , and the sum of the squares of the electric 

P 2 0 2 2 

and quadrupole form factors of the deuteron, g£ + — p Gp. . 

d 9 y d 

Because of the smallness of n in our range of momentum 

transfers (n = 2.2 x 10 - ^ at q^ = 0.80 F~^) the quadru- 
max 

-4 

pole term contributes only about one part in 10 to the 

form factor and can be safely neglected. 

The end result of the experiment is then the ratio 
2 2 

Gg /Gg from which the ratio Gg /Gg is extracted with the 
d p n p 

help of equation (2-128). 

In summary, the quantity evaluated from the data is 
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sc 
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sc 



x n ? x cP x °Mott 



x K P x kJ 
s b 



(4-5) 
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where C represents the magnetic correction, °]y[ 0 tt ^he Mott 
cross section, and the indices p,d refer to proton or 
deuteron, respectively. 

C . EXPERIMENTAL ERRORS 

Errors encountered in the measurements are of both 
statistical and systematic nature. Systematic errors, such 
as the value of the solid angle, counting rate correction, 
target angle, scattering angle, and SEM efficiency cancelled 
out since they were measured under the same experimental 
conditions. The absolute error in the energy of the incom- 
ing electrons is about 0 . 2 %. Since the Mott cross section 
for the deuteron and the proton are evaluated at nearly the 
same energy, this error becomes insignificant. The same is 
true for the errors in the radiative corrections if they are 
evaluated over roughly the same energy cut off AE (Equation 

4 - 3 ). 

The only significant error of a systematic nature is 
the error in the target thickness and density, pt. Because 
of the thinness of the targets, it could be established to 
only 0 . 5 %. 

The statistical error consists of two parts, the actual 
counting error and the error in the efficiency correction. 

It is in principle possible to keep these errors arbitrarily 
small by taking a very high number of counts. Restrictions 
in counting speed as posed by the electronic equipment and 
time limitation in the use of the accelerator pose a 
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practical limit. All but two peaks contain ^10,000 counts, 

which would make their statistical error to be 0.5$. The 

required carbon subtraction introduces an additional error, 

depending on the number of carbon counts under the peak area. 

The error in the efficiency corrections is quite ini- 

formly 0.5$ for the Naval Postgraduate School runs, where 

roughly the same number of counts have been taken to 

establish the relative efficiencies. In the Stanford runs 

the error varies from 0.27$ to 0.65$. 

2 2 

For each measurement of G^. /G^ , the error is then the 

E d E p 

square root of the quadratic sum of the counting, efficiency, 
and target error. 

One half of this percentage error enters the ratio 

G F /Gp and, therefore, the calculation of G„ /G„ . This 
a p n p 

is the error quoted in Table VII. 
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V. COMPARISON TO THEORY AND DISCUSSION 



A. EXPERIMENTAL RESULTS 

The experimental data are presented in Table VI. The 

error quoted consists of counting errors (including the 

carbon subtractions and efficiency corrections) and an error 

of 0.5% in the target thickness. 

Table VII shows four sets of values of Gg /Gg as 

n p 

obtained by: (a) the Feshbach-Lomon wave function with 

relativistic corrections; (b) the Feshbach-Lomon wave func- 
tion without relativistic corrections; (c) the Partovi 
wave function with relativistic corrections, and; (d) the 
Partovi wave function without relativistic corrections. 

Values for the relativistic correction AGg and, as an 

n 

example, the structure factor F^ for the Feshbach-Lomon 
wave function are also given. In each of the columns are 

p 

shown the slopes ( dGg /Gg )/dq to illustrate the rapid 

n p 

decrease in slope when going from a relativistically 
corrected Feshbach-Lomon wave function to an uncorrected 
Partovi wave function. 

The b' fit of de Vries [28] for the absolute proton 

charge form factor was used to extract Gg from the Gg /Gg 

n n p 

ratio. This fit is in very good agreement with the absolute 

measurements of Gg by Drickey and Hand [1] . Table VIII 

P 

shows the results of the b' fit for the two extreme cases 
presented in Table VII. The slopes as determined by the 
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the two sets of data are also given for comparison with 
the neutron-electron interaction slope. 

B. DISCUSSION 

The Figures 19 and 20 illustrate the results tabulated 
in Tables VII and VIII respectively. It can be seen from 
the slope of Figure 20 that the Partovi wave function with- 
out relativistic corrections does not significantly differ 
from the data of Drickey and Hand pictured in Figure 10(a). 
That is, they both essentially agree that the slope is very 
small or zero. On the other hand the slope for the Feshbach- 
Lomon wave function with relativistic corrections approaches 
closely the value of the slope for the neutron-electron 
interaction . 

Barring unexpected fluctuations in the proton form 
factor G„ in the range of momentum transfers of interest 

hi 

P 

here it can be concluded that: 

(a) The Feshbach-Lomon wave function together with the 

relativistic corrections AGg , removes the apparent discrep- 

n 2 

ancy between the neutron-electron slope at q =0 and the 

slope given from values of G p (obtained by electron scat- 

n ■ -2 

tering) in the range 0 . 10 <_ q^ <_ 0 . 80 F 

(b) Even within the relatively large errors propagated 

into Gg , the Partovi wave function with relativistic 
n 

corrections applied, is in disagreement with the neutron- 
electron interaction slope. 



112 



(c) With the present small experimental uncertainties. 
It is reasonable to state that this experiment should be 
added to the group of experiments which must be explained 
by any deuteron theory. That is, a proper model for the 
deuteron must be found that will predict the slope as 
determined by the neutron-electron scattering experiments. 

The important qualitative conclusion from these experi 
mental results is that the neutron has a non-zero charge 
form factor. This implies that there is a charge distribu 
tion within the neutron. The neutron thus must have a 
layered charge structure with the outermost part of the 
charge distribution being negative in sign. 
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TABLE VIII 



SLOPE OP THE CHARGE FORM FACTOR 

dGg /dq 2 
n 



q 2 

(F -2 ) 


G t? a 

e p 

de Vries 
b’ fit 


p FL 

g e 

n 


g e +ag e 

n n 


r P 

g e 

n 


0.10 


0.9881 


-0.0027 


-0.0020 ± 0.0054 


-0.0038 ± 0.0054 


0,20 


0.9765 


+0.0039 


+0.0051 ± 0.0031 


+ 0.0017 ± 0.0031 


0.30 


0.9652 


-0.0006 


+ 0.0012 ± 0.0035 


-0.0037 ± 0.0035 


0, 4o 


0.9540 


+0.0077 


+ 0.0103 ± 0.0036 


+ 0.0036 ± 0.0036 


0.50 


0.9*131 


+0.0037 


+ 0.0067 ± 0.0074 


-0.0015 ± 0.0074 


0,60 


0.9321 


+0.0011 


+ 0,0047 ± 0.0066 


-0.0042 ± 0.0066 


0.80 


0.9H5 


+0.0138 


+ 0.0184 ± 0.0059 


+0.0076 ± 0,0059 
















Slope : 


0.0179 ± 0.0036 


0.0026 ± 0.0036 



Slope of n-e interaction = 0.0193 - 0.0004 



a: No error has been applied to Gg . 
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Figure 19 
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Figure 20 



APPENDIX A 



THE NEUTRON CHARGE FORM FACTOR FROM THE 
NEUTRON-ELECTRON INTERACTION 



In the n-e experiments what is measured is a scattering 

amplitude of a neutron by a bound electron. In the Born 

approximation for neutrons of very long wavelength this is 

related to the volume integral of V(r), the interaction 

potential [253 . The magnitude of the n-e interaction is 

expressed in the form of an equivalent potential Vq over a 

e 2 

sphere of radius r n = the classical electron radius, 

u m c d 

G 

so that the volume integral of Vq is the experimental value 
of the potential. The potential Vq is simply a convention 
since the classical electron radius plays no fundamental 
role in the problem. Thus, 



/■ 



V(r)dr 



V(r ) = V 



0 



const . 




r 3 
3 0 



V, 



(A-l) 



Using Poisson's equation and the definition, of r.m.s. radius 
p 2 ->■ 

<r > 5 /r p(r)dr, (A-l) becomes 

Jv( r)dr = Jr 2 p(r)dr = <r 2 > (A-2) 

Equating the right hand sides of (A-l) and (A-2) yields, in 
rationalized MKS units 



< r 2 > = r V (eV) 1.6 x 10 -19 

3 X H TT£ 3 o o 



. x2 



where e is the permitivity of free space and in (Fermi) 
o 
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( A— 3 ) 



0 V eV x r x 8ire 0 
^ 2 . o o o „2 

< r > = — p 



e x 10 



-30 



<r 2 > = 3.1^ x 10 5 V q (eV) F 2 . 

From (2-90) the expression for the charge form factor is 
given by 



2 . 2 . 

1 ■■ 

n 



( A-4 ) 



so that the slope at the limit q 0 gives 



'dG T 



n 



dq‘ 



-<r 2 > 



2 

q ^->0 



-5-23 xlO" 6 V o (eV) F 2 (A-5) 



Krohn and Ringo's [16] latest value is V o = -3720 eV, hence. 



dG 



E 



n 



dq' 



q 2 +0 



= +0.0193 F‘ 
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APPENDIX B 



PROPERTIES OF THE TENSOR OPERATOR S 12 
Several experimental facts, like the nonadditivity of 



of the quadrupole moment in the deuteron, led to the 
suggestion that the nature of the nuclear force is 
partially noncentral, that is, it contains a partly "tensor" 
force . 

Require that the tensor operator be symmetric in 
particles 1 and 2, and that it be invariant under the parity 
transformation (space reversal). Let r = be the 

relative coordinates of particle 1 and 2. The particles 
have spins given by and o ^ where a's are the Pauli spin 
matrices (an explicit representation) 



the magnetic moments (i.e., + P n ^ P^) and the presence 




(B-l ) 



These matrices obey Pauli matrix algebra 




26 . i,J = x,y,z 



( B-2 ) 



a x a 



2ia 



( B— 3 ) 




(B-*») 



(a • it) (a •§) = ()t*$) + ia*(^x§) 

Application of the parity operator P to r and c yields 



(B-5) 



Pr -r. 



Pa -* a . 
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Thus the vector r can appear only an even number of times 
in the tensor operator. Further, since the tensor operator 
is symmetric in 1 and 2, and a 2 should appear in some 
bilinear combination. With (B-2,3,^j5) in mind it is found 
that the only linearly independent quantities which can be 
made from r, and o 2 are 

( o 1 • ct 2 ) and (a 1 *r)( 0 2 *r) (B-6) 



The traditional definition then for the tensor operator 
S^ 2 is 

S 12 = 3(o’ 1 *r) (o 2 *r) - (a-^c^) (B-7) 



where r = r/ | r | . 

The tensor operator (B-7) has some interesting proper- 
ties, namely <S^ 2 > over all directions vanishes. Consider 



< ( a.. *r) (a. 




sin 6- d9 d<J> 



(B-8 ) 



Take along z-axis , o 2 in the x-z plane making angle a 
with the z axis, and r a vector making angles ©-, and <J> . 
Then 



°1 = °l z 

-Jy. 

°2 ~ a 2^ s ^- nax + cosaz ) 



A 

r 



A A A 

sin 9- cos <J> x + sin 9- sin <j> y + cos 9- z 



(B-9 ) 
(B-10) 
(B-ll ) 



Putting (B-9, 19, 11) in (B-8) and integrating over 9- and <J> 
yields 
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(B-12 ) 



< (c^ *r) (o 2 • r)> = ^ cos a = (c^a^ 



Hence it follows that <S^ 2 > = 0. 
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APPENDIX C 



CALCULATION OP DEUTERON MAGNETIC 
AND QUADRUPOLE MOMENTS 



The neutron and proton magnetic moments y ^ , and y^ have 

been quite accurately determined } as has the deuteron' s 

magnetic moment y^. The difference (y + y^) - y^ = 0.0222nm 

is small but not zero and much larger than the experimental 

error. A possible explanation is that the deuteron spends 

a small part of the time in some state other than the 

q 

symmetric state. The task is to find the correct mixture 

of states that can account for the deuteron' s moment. 

S 3 

Recall the angular wave functions for the and D^ 

states 



'rVoii- 



v° 1 

= Y o X 1 



(2-56) 



, D * 
U l’ 



211 



- [jL -1 _ /_3_ y 1 + I— y° v"*- (2-57) 

1 V 5 X 2 X 1 \10 2 X 1 \10 2 X 1 ^ 0 



Using (2-5 6 ) and (2-57) construct the deuteron wave function 
Thus 



3s i : *011 (?) 
3 d 1 : * 211 ^ 





(C-l) 
( C— 2 ) 



The nuclear magnetic 



z 




k=l 



moment operator can be written 

A 

+ ®p^zk^ + ®n ^ ^zk 
k=z+l 



(C-3) 



12 ^ 



t h 

where is the orbital angular momentum of the k proton, 

X_ 

the spin angular momentum of the k proton (or neutron) 
and where g and g are the gyromagnetic ratios of the 

ir * 

proton and neutron, g^ = 5.59, g R = -3.83. (The carot over 
y_ denotes y_ as an operator.) In the c-m system and for 
the deuteron ( C— 3 ) becomes 



y = rr i + g S + g S 
z 2 z b p zp & n zn 



(C-M 



where the orbital angular momentum associated with the 
proton is half of the relative orbital angular momentum. 

The magnetic moment is found by taking the expectation value 



of (C-M between the state (f> 



m 



£SJ 



'ft 



„ v. , , >#m * ,m ,.m i ^ i , m 

<y z > = dr d) £SJ y z ^ S J = K hsj I y z ' *£SJ > 



(C-5) 



The following formula are useful in evaluating (C-5) 
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zp X 0 2 X 1 



^zn X 1 



1 x° 

2 X 1 
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Hence 




The orthonormality of the spherical harmonics Y integrated 

over the solid angle, and the orthonormality of the spin- 

m„ 

wave functions x c applied to (C-7) yields 

o 



3 





g o 6 n T 2 • S n + S n 

- / u^(r)dr = -£ 2. 

2 J 2 



(C-8) 



= 0.8797 nm 



oo 

S 1 <4> 211 1 ’"z I ^211 >= [f - T (s p +g n ) ]/ w2(r,dr = 



(C-9) 



[f-J (« P + B n ) ] = 0 ' 3101 nm 



In evaluating the radial integrals the assumption was 
made that the deuteron was in a definite state, either S or 
D but not both. But the deuteron magnetic moment is 
0.8573 nm, that is, neither of the two values above. Thus 
consider a mixture of the two states. 

The total wave function for a combination of states is 



<p = <P + 
' s 



♦d = ♦oil + 



211 



(C-10) 
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Normalization gives 



<<f>|(f» = P g + P D = 1 (C-ll) 

00 2 °° 2 

where P = / u (r)dr, P n = / w (r)dr are the individual 
S 0 D 0 

probabilities that the deuteron is in an S or D state. Thus 

U 3 s +3 n = ! <g p +S n )P s + [|-|(Sp +g n , ] P D (C - 12 > 

1 1 

Equations (11) and (12) are two simultaneous equations in 
the two unknowns P^ and Pp. Using the experimental value 
of y^ for y in (C-12) the solution for P g and P D are 

P D = *i%, P s = 96% 

The quadrupole moment can be found in the same manner 
as the magnetic moment, 

Q d = - r = (3 cos^O- - l)r (C-13) 

In the c-m system, the distance of the neutron and proton 
from the center of mass is r/2 , and only the proton con- 
tributes to Q^, Hence, 

< Q d > = ^-^>*(3 cos^U- - l)4>r^dr (C-l4) 

In terms of spherical harmonics this becomes. 
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The first term represents the pure “'S^ contribution and when 



integrated over the solid angle it is zero because of 
spherical symmetry. The next terms in parentheses are 
hermitian conjugates of each other and the total can be 
written as 



yfW f uw1 J on Y 2 V211 d? 



which after angular integration gives 



1 

/50 



a 

/ 



uwr dr. 



The third term comes from the state and gives 



- -f 

20 J 



2 2, 
w r dr . 



Hence 






<x 

10 J 



2 , 1 
uwr dr - — 



oc 

/ 



2 2 , 

w r dr 



( C-15 ) 



To go further, explicit solutions for the radial wave func- 
tions u(r) and w(r) are needed. Some examples of for 
specific numerical solutions u and w are listed in Table III. 
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APPENDIX D 



THE BREIT OR BRICKWALL FRAME 
The Breit frame is a Lorentz frame in which either 

4- 4 | 4- 4- f 

k-j_ + k ^ or, p^ + vanish, where in an elastic two-body 

4" 4" 4* t 

collision, k^, are incident particle 3-momentum and k 2 , 
P 2 are outgoing particle 3-momentum. The somewhat odd 
notation stems from traditional techniques in elastic 
scattering theory and dispersion relations, particularly 
the Mandelstam representation. Consider Figure D-l. The 
convention is that all four-momentum are ingoing. This has 
the advantage that one may consider any two of the four to 
be the incoming particles and the other two as outgoing 
ones; the physical momentum of the outgoing particles is 
then the negative of the one used in Figure D-l. 




The physical momentum of an outgoing particle is denoted 
by a prime. For example, if the particle corresponding to 
the arrow with k 2 is outgoing, it has physical momentum 
k 2 = -k2* 
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Apply a Lorentz transformation to the scattering event 
such that It + = 0. Thus, four-momentums and have 

the form (h = c = 1) 

k 1 = (w,k), k' 2 = ( w , -k ) ( D-l ) 

From conservation of energy it follows that the energies of 
the p particle before and after the collision, and e 2 
must be equal. Hence | | = | p ^ I and 

P x = P 2 = ( c >v' 2 ) (D-2) 



with 




Since scattering is elastic k^ - k^ = P 2 - p-^ and it follows 

k 2 - k 2 = (0,2k) = p 2 - Pi = ( 0 ,Pg ~ P 2 ) (D-3) 

2k the three momentum transfer. Equation (1) (2) and (3) 
yield the following picture. Both particles seem to be 
reflected on a hard wall, the k particle perpendicularly and 
the particle at some angle 




Figure D-2 
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In the c-m system the energies e and w of the p and k 
particles are not independent, but they are in this system. 
Therefore they are convenient to use as variables. In 
this system the four momentum transfer is very simple: 

q 2 = (k- L - ) 2 = (0,21c) 2 = - 1 21c | 2 =-q 2 

i.e., it is the three momentum transfer. 
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